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This  DEPSCoR/AFOSR  sponsored  research  grant  has  been  used  to 
continue  the  long-term  support  of  our  graduate  level  research  program 
on  the  asymptotic  description  of  ultrawideband  signal,  ultrashort  pulse 
electromagnetic  wave  propagation  in  causally  dispersive  media  and 
waveguiding  systems,  extending  it  to  include  the  properties  of  spatial 
inhomogeneity.  Our  long-term  research  goal  in  this  important  area  is  to 
develop  a  rigorous,  uniform  asymptotic  description  of  ultrawideband 
electromagnetic  pulsed  beam  propagation,  reflection,  and  transmission 
phenomena  in  both  lossy  dielectric  and  conducting  dispersive  media  that 
may  also  exhibit  additional  complicating  behavior  such  as  either  spatial 
inhomogeneity  or  spatial  dispersion  that  is  valid  for  arbitrarily  short  rise¬ 
time  pulses.  The  results  of  this  analysis  have  direct,  meaningful 
application  to  the  analysis  and  design  of  low-observable  surfaces  (for 
stealth  airframes)  and  ultrawideband  radar  systems  (for  observing 
stealth  airframes),  the  remote  detection  of  buried  structures  (such  as 
landmines  and  IED’s),  ionospheric  pulse  propagation  (for  remote  sensing 
from  an  orbiting  satellite),  as  well  as  the  problem  of  ultrawideband 
electromagnetic  pulse  exposure  of  biological  tissues.  Of  further  interest  is 
the  application  of  this  theory  to  undersea  communications  systems, 
terahertz  optical  communication  and  integrated  optics  systems,  and  the 
remote  sensing  of  geophysical  structures. 

This  AFOSR  sponsored  research  program  has  served  to  provide  partial 
support  for  my  research  colleague  Dr.  Natalie  Cartwright  as  a  Research 
Associate  Professor  in  the  Department  of  Mathematics  at  the  University 
of  Vermont  during  the  final  two  years  of  this  grant  (2006-2007)  following 

20071101407 


the  successful  completion  of  her  Ph.D.  degree  in  2004.  A  copy  of  her 
dissertation  "Uniform  Asymptotic  Description  of  the  Unit  Step  Function 
Modulated  Sinusoidal  Signal"  was  reprinted  in  the  University  of  Vermont 
ECE/04/05-01  Research  Report  that  was  submitted  with  my  2004  Final 
Technical  report.  The  paper  “Uniform  Asymptotics  Applied  to 
Ultrawideband  Pulse  Propagation”,  based  upon  her  dissertation,  is 
scheduled  to  be  published  in  SIAM  Review  in  November  of  this  year.  This 
important  work  completes  the  uniform  asymptotic  description  of 
dispersive  pulse  propagation  in  a  Lorentz  model  dielectric,  a  problem  that 
was  initiated  by  Arnold  Sommerfeld  in  1907  and  by  Leon  Brillouin  in 
1914,  was  revisited  by  Baerwald  in  1930,  Stratton  in  1941,  Oughstun 
and  Sherman  in  1988,  and  finally  completed  by  Cartwright  and 
Oughstun  in  2007,  a  mere  century  after  Sommerfeld's  first  publication 
regarding  the  mathematically  proper  solution  to  this  fundamental 
problem.  The  asymptotic  solution  to  the  ionospheric  propagation 
problem,  including  the  effects  of  spatial  inhomogeneity,  is  currently  being 
pursued  by  Professor  Cartwright  at  SUNY-New  Paltz  under  a  separate 
AFOSR  Young  Investigator  Research  award.  This  challenging  research 
direction  is  a  result  of  the  research  developed  under  this  grant  and  is 
being  conducted  with  my  collaborative  support. 

As  a  product  of  this  funded  research,  the  following  papers  were 
published: 

1 .  K.  E.  Oughstun  and  N.  A.  Cartwright,  “Physical  Significance  of  the 
Group  Velocity  in  Dispersive,  Ultrashort  Gaussian  Pulse  Dynamics," 
Journal  of  Modem  Optics  52,  8,  1089-1 104  (2005). 

2.  K.  E.  Oughstun,  “Dynamical  Evolution  of  the  Brillouin  Precursor  in 
Rocard-Powles-Debye  Model  Dielectrics,"  IEEE  Transactions  on 
Antennas  and  Propagation  53,  5,  1582-1590  (2005). 

3 .  K.  E.  Oughstun,  "Several  Controversial  Topics  in  Contemporary 
Optics:  Dispersive  Pulse  Dynamics  and  the  Question  of  Superluminal 
Pulse  Velocities,"  in  Masters  in  Optics:  A  Tribute  to  Emil  Wolf  (SPIE 
Press,  2005)  pp.42 1-453. 

4.  K.  E.  Oughstun  and  R.  A.  Albanese,  "Magnetic  Field  Contribution  to 
the  Lorentz  Model,"  Journal  of  the  Optical  Society  of  America  A  23.  7 
1751-1756  (2006). 

5.  N.  A.  Cartwright  and  K.  E.  Oughstun,  “Uniform  Asymptotics  Applied 
to  Ultrawideband  Pulse  Propagation,”  SIAM  Review  49,  4,  xxx-xxx 
(2007). 

6. N.  A.  Cartwright  and  K.  E.  Oughstun,  “Ultrawideband  Pulse 
Propagation  through  a  Lossy  Plasma,”  Radio  Science  (in  preparation). 

Copies  of  these  published  papers  are  attached  at  the  back  of  this  report. 
Reprints  have  already  been  sent  to  both  Dr.  Arje  Nachman  at  AFOSR  and 
to  Dr.  Richard  Albanese  at  Brooks  City  Base  at  the  time  each  publication 
appeared. 
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In  addition,  the  following  conference  and  workshop  presentations  were 

given  during  the  period  of  this  grant  support: 

1 .  K.  E.  Oughstun,  R.  A.  Albanese  and  J.  Penn,  “Trapezoidal  Envelope 
Pulse  Dynamics  in  Debye  Model  Dielectrics,"  2005  AFOSR 
Electromagnetics  Workshop  (January  5-7,  2005,  San  Antonio,  Texas). 

2 .  N.  A.  Cartwright  and  K.  E.  Oughstun,  “Uniform  Asymptotic 
Description  of  the  Main  Signal  Arrival  in  Dispersive  Pulse 
Propagation,”  2006  AFOSR  Electromagnetics  Workshop  (January  11- 
13,  2006,  San  Antonio,  Texas). 

3 .  N.  A.  Cartwright  and  K.  E.  Oughstun,  “Uniform  Asymptotic 
Description  of  the  Main  Signal,"  2006  USNC/URSI  National  Radio 
Science  Meeting,  (January  4-7,  2006,  University  of  Colorado,  Boulder, 
CO). 

4 .  N.  A.  Cartwright  and  K.  E.  Oughstun,  “Ultrawideband  Pulse 
Propagation:  The  Signal  Contribution,"  AMS  Joint  Mathematics 
Meeting  (San  Antonio,  Texas). 

5.  N.  A.  Cartwright  and  K.  E.  Oughstun,  “Uniform  Signal  Contribution  of 
the  Step  Function  Modulated  Sine  Wave,"  Progress  in 
Electromagnetics  Research  Symposium  (PIERS  2006)  (March  26-29, 
2006,  Cambridge,  MA). 

6 .  N.  A.  Cartwright  and  K.  E.  Oughstun,  “Uniform  Asymptotic 
Description  of  the  Main  Signal  Contribution  in  Dispersive  Pulse 
Propagation,"  Northern  Optics  2006  (June  14-16,  2006,  Bergen, 
Norway). 

7.  K.  E.  Oughstun  and  N.  A.  Cartwright,  “Dynamical  Pulse  Evolution  in 
the  weak  Dispersion ,  Limit,"  Eightennth  Annual  AFOSR 
Electromagnetics  Workshop  (January  9-11,  2007,  San  Antonio, 
Texas). 

8.  K.  E.  Oughstun,  “Propagation  through  Dispersive,  Absorptive  Media,” 
AFOSR  Workshop  on  Electromagnetic  Wave  Propagation  through 
Challenging  Media  (May  8-9,  2007,  Hanscom  AFB). 

9 .  N.  A.  Cartwright  and  K.  E.  Oughstun,  “Ultrawideband  Pulse 
Penetration  in  a  Debye  Medium  with  Static  Conductivity,”  Fourth 
IASTED  International  Conference  on  Antennas,  Radar,  and 
Propagation  (May  30- June  1,  2007,  Montreal,  Canada). 

10.  N.  A.  Cartwright  and  K.  E.  Oughstun,  “Ultrawideband  Pulse 
Penetration  in  an  Isotropic  Collisionless  Plasma,"  2007  CNC/USNC 
North  American  Radio  Science  Meeting  (July  22-26,  2007,  Ottawa, 
Canada). 

11.  K.  E.  Oughstun,  “Optimal  Pulse  Penetration  through  Dielectric 
Barriers,”  2007  CNC/USNC  North  American  Radio  Science  Meeting 
(July  22-26,  2007,  Ottawa,  Canada). 

12.  K.  E.  Oughstun  and  N.  A.  Cartwright,  “Ultrashort  Electromagnetic 
Pulse  Dynamics  in  the  Singular  and  Weak  Dispersion  Limits," 
Progress  in  Electromagnetics  Research  Symposium  (PIERS  2007) 
(August  27-30,  2007,  Prague,  Czech  Republic). 
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Finally,  the  first  volume  of  my  two  volume  work  Electromagnetic  and 
Optical  Pulse  Propagation  1:  Spectral  Representations  in  Temporally 
Dispersive  Media  on  time-domain  electromagnetics  was  published  in  the 
Springer  Series  in  Optical  Sciences  in  2006.  The  Preface  in  this  book 
states  that  "the  critical,  long-term  support  of  this  research  by  Dr.  Arje 
Nachmann  at  the  Physics  and  Electronics  Directorate  of  the  United 
States  Air  Force  Office  of  Scientific  Research  is  gratefully  acknowledged.” 
The  second  volume  Electromagnetic  and  Optical  Pulse  Propagation  2: 
Temporal  Pulse  Dynamics  in  Dispersive,  Attenuative  Media  is  near 
completion. 
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Physical  significance  of  the  group  velocity  in  dispersive,  ultrashort 

gaussian  pulse  dynamics 

KURT  E.  OUGHSTUN*  and  NATALIE  A.  CARTWRIGHT 

College  of  Engineering  and  Mathematics, 

University  of  Vermont,  Burlington,  VT  05405,  USA 

(Received  16  February  2004;  in  final  form  26  April  2004) 


The  properties  of  ultrashon  gaussian  pulse  propagation  in  a  dispersive,  atienua- 
live  medium  are  reviewed  with  emphasis  on  the  pulse  velocity.  Of  particular 
interest  is  the  group  velocity  whose  physical  interpretation  loses  meaning  in 
causally  dispersive  materials  as  the  temporal  pulse  width  decreases  into  the 
ultrashort  pulse  regime.  A  generalized  definition  of  the  group  velocity  that  applies 
to  ultrashort  pulses  in  causally  dispersive  materials  is  provided  by  the  centroid 
velocity  of  the  pulse  Poynting  vector  whose  properties  are  described  here. 
In  particular,  it  is  shown  that  this  physical  velocity  measure  approaches  the 
group  velocity  for  any  value  of  the  initial  pulse  carrier  frequency  and  at  any  fixed 
value  of  the  propagation  distance  in  the  limit  as  the  initial  pulse  width  increases 
indefinitely.  This  then  provides  a  convenient  measure  for  determining  when  the 
group  velocity  approximation  is  valid. 


1,  Introduction 

The  interrelated  problems  of  dispersive  wave  propagation  and  the  group  velocity 
have  a  long  and  involved  history.  In  1839,  Sir  William  R.  Hamilton  [1]  considered 
dispersive  wave  propagation  as  a  coherent  superposition  of  monochromatic  scalar 
wave  disturbances,  introducing  the  concept  of  group  velocity  and  comparing  it  to  the 
phase  velocity.  Lord  Rayleigh  (2]  mistakenly  attributed  the  original  definition  of  the 
group  velocity  to  Stokes  [3]s  stating  that  ‘ when  a  group  of  waves  advances  into  mill 
water ,  the  velocity  of  the  group  is  less  than  that  of  the  individual  waves  of  which  it  is 
composed:  the  wavejr  appear  to  advance  through  the  group,  dying  away  as  l hey 
approach  its  anterior  limit.  This  phenomenon  was,  I  believe ,  first  explained  by  Stokes, 
who  regarded  the  group  as  formed  by  the  superposition  of  two  infinite  trains  of  n-avej, 
of  equal  amplitudes  and  of  nearly  equal  wavelengths §  advancing  in  the  same  direction  . 
Rayleigh  [4]  then  used  these  results  to  explain  the  difference  in  observability  of  the 
phase  and  group  velocities. 

The  distinction  between  signal  and  group  velocities  originated  in  the  early 
research  of  Voigt  [5,  6]  and  Ehrenfest  [7]  on  elementary  dispersive  waves,  and  by 
Laue  [8]  who  first  considered  dispersive  wave  propagation  in  a  region  of  anomalous 
dispersion  where  the  absorption  is  both  large  and  strongly  frequency  dependent.  The 
distinction  between  front  and  signal  velocities  was  then  considered  by  Sommerfeld 
J9,  10]  who  proved  that  no  signal  could  travel  faster  than  the  vacuum  speed  oflight 
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The  first  experimental  measurement  of  the  signal  velocity  was  attempted  by  Shiren  [2 1  ] 
in  1962  using  pulsed  microwave  ultrasonic  waves  within  a  narrow  absorption  band, 
A  more  detailed  analysis  of  these  experimental  results  by  Weber  and  Trizna  [22] 
indicated  that  the  velocity  measured  by  Shiren  was  in  reality  that  for  the  first 
precursor  and  not  the  signal.  Subsequent  research  by  Handetsman  and  Blcistein  [23] 
in  1969  provided  a  uniform  asymptotic  description  of  the  arrival  and  initial  evolution 
of  the  signal  from.  The  first  experimental  measurements  of  the  precursor  fields 
originally  described  by  Sommerfeld  [10]  and  Bnlloum  [12,  13]  were  then  published 
by  Pleshko  and  Palocz  [24]  who  first  referred  to  the  first  and  second  precursors  as 
the  Sommerfeld  and  Brillouin  precursors,  respectively.  Their  experimental  results 
established  the  physical  propriety  of  the  asymptotic  approach. 

The  establishment  of  the  equivalence  between  the  group  velocity  and  the  energy 
transport  velocity  in  loss-free  media  [25-27]  provided  a  physical  basis  for  the  group 
velocity  in  lossless  systems.  The  precise  formulation  of  the  quasimonochromatic 
or  slowly-varying  envelope  approximation  by  Born  and  Wolf  [28]  in  the  context  of 
partial  coherence  theory  then  completed  the  mathematical  and  physical  basis  for  the 
group  velocity  approximation,  which  was  then  generalized  [29]  and  extended  [30,  31] 
to  any  order  of  dispersion.  The  quasimonochromatk  or  slowly-varying  envelope 
approximation  is  a  hybrid  time  and  frequency  domain  representation  [32]  in  which 
the  temporal  field  behaviour  is  separated  into  the  product  of  a  temporally  slowly 
varying  envelope  function  and  an  exponential  phase  term  whose  angular  frequency  is 
centred  about  some  characteristic  frequency  of  the  pulse.  The  envelope  function  is 
assumed  to  be  slowly  varying  on  She  time  scale  A tc  **  \/a>c,  w'hich  is  equivalent  [33] 
to  the  assumption  that  its  spectral  bandwidth  Aw  satisfies  the  inequality  Aw/wc  <£  I. 
The  frequency  dependence  of  the  wavenumber  may  then  be  approximated  by  the  first 
few  terms  of  its  Taylor  series  expansion  about  the  characteristic  pulse  frequency 

with  the  assumption  [30^32]  that  improved  accuracy  can  always  be  obtained 
through  the  inclusion  of  higher-order  terms;  this  assumption  has  been  proven 
incorrect  [34-35],  optimal  results  being  obtained  using  either  the  quadratic  or  the 
cubic  dispersion  approximation  of  the  wavenumber. 

The  description  of  the  velocity  of  energy  transport  through  a  causally  dispersive 
medium  [12,  13]  was  reinvestigated  by  Schulz-DuBois  [36]  in  1969  and  finally  by 
Loudon  [37]  in  1970  who  provided  a  correct  description  of  the  energy  velocity 
in  a  single  resonance  Loreni2  model  dielectric.  This  description  and  its  extension  [38] 
to  a  multiple  resonance  Lorentz  model  dielectric  showed  that  the  energy  velocity  and 
group  velocity  are  different  in  the  region  of  anomalous  dispersion.  Based  upon  this 
critical  result,  Sherman  and  Oughstun  [39,  40]  then  presented  a  complete  physical 
description  of  dispersive  pulse  dynamics  in  a  causally  dispersive  medium  in  terms  of 
the  energy  velocity  and  attenuation  of  time-harmonic  waves.  This  description 
reduces  to  the  approximate  group  velocity  description  in  the  limit  as  the  material 
loss  goes  to  zero,  A  precise  description  of  the  signal  velocity  was  then  given  by 
Oughstun  and  Sherman  [41]  in  1988  in  connection  with  the  modern  asymptotic 
theory  of  dispersive  pulse  propagation  [41-44],  The  physical  propriety  and  related 
observability  of  this  signal  velocity  definition  was  then  demonstrated  (44,  45]  through 
a  numerical  experiment,  thereby  completing  the  physical  interpretation  of  the 
asymptotic  description. 

It  is  clear  that  a  more  physically  meaningful  pulse  velocity  measure  needs  to  be 
considered  in  order  to  accurately  describe  the  complicated  pulse  evolution  that 
occurs  in  ultrashort  dispersive  pulse  dynamics.  One  possible  measure  is  given  by  the 
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In  the  asymptotic  theory  of  dispersive  pulse  propagation  [12,  44]  the  integral 
representation  given  in  equation  (1)  is  expressed  as 

a(2,  /)  =  ^  \m  exp  |(Az/c)0{w,6i)]  d  w  (5) 

with  complex  phase  function 

$(&>,  Q)  =  i  |^(cu)Ar  -  iu/J  ~  icu{«(iu)  —  0\  (6) 

and  non-dimensional  space-time  parameter  $  =  erf  Az.  The  fact  that  this  exact 
integral  representation  of  the  propagated  optical  wave  field  satisfies  relativistic 
causality  is  expressed  by  the  following  theorem  (originally  proved  by  Sommerfeld 
[9,  10]  for  a  Heaviside  unit  step  function  modulated  signal  in  a  single  resonance 
Lorentz  model  dielectric  and  later  extended  [39,  44]  to  an  arbitrary  plane  wave  pulse 
in  a  general  causally  dispersive  medium): 

Sommerfetds  Theorem:  If  /(f)  —  0  for  all  t  <0  and  if  JH  {ia>[rc(cu)  -  0]}  — ►  — oo  as 
\to\  — ►  oo  with  3  {co\  =  \a>\  sin(tf)  for  arbitrarily  small  >  0  for  all  9<  h  then 
A(z,  t)  =  0  for  all  Az  >  G  when  0<  L 

This  precise  statement  of  the  luminal  arrival  of  the  signal  front  then  proves  that  any 
information  that  may  be  present  in  the  signal  will  follow  at  some  later  space- time 
point  with  0  >  I . 


3,  The  velocities  of  a  gauss!  an  envelope  pulse  in  Lorentz  mode)  dielectrics 

The  group  velocity  approximation  is  a  hybrid  time-frequency  domain  representation 
[32]  in  which  the  pulse  is  separated  into  the  product  of  a  slowly-varying  envelope 
and  an  exponential  phase  with  angular  frequency  centred  about  a  characteristic 
frequency  of  the  initial  pulse.  With  the  initial  pulse  at  the  plane  z  =  r0  given 
by  /(f)  =  i/(/)  sin  (a»c/  +  \j/)  with  envelope  i/(f)  and  fixed  carrier  frequency 
the  propagated  plane  wave  pulse  is  then  given  by  the  Fourier-Laplace  integral 
representation  [44] 

A[z,  f)  —  2“  SR  |i  exp(-i^)  J  u(co  -  (dc)  exp  jij^:(uj)Ar  -  eu/J  J  dci>J  (7) 

for  all  Ar  >  0,  where  yp  —  0,  jt/2  for  either  a  cosine  or  sine  wave  carrier,  respectively. 
Here  u{a>)  is  the  temporal  angular  frequency  spectrum  of  the  initial  pulse  envelope 
function  w(f).  Because  of  its  central  importance  in  optics,  a  unit  amplitude  gaussian 
envelope  pulse  is  considered  here,  where 

u(t)  =  zxp(-r/T2),  (8) 

with  initial  full  pulse  width  27  >  0  measured  at  the  exp{—  1)  amplitude  points. 
In  a  multiple  resonance  Lorentz  model  [11]  dielectric  the  complex  index  of  refraction 
is  described  by  [44] 
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In  particular,  the  cubic  dispersion  approximation 

m  +  W«  ~^)+~ -  *cf  +  -  ft*)3  ( 1 4) 

is  also  used  because  the  cubic  term  introduces  a  small  degree  of  asymmetry  imo  the 
propagated  pulse.  However,  recently  published  research  [34,  35]  has  established  that 
optimal  results  in  the  global  sense  are  obtained  either  with  the  quadratic  dispersion 
approximation  (10)  or  with  the  cubic  dispersion  approximation  (14). 

The  phase  velocity  vp(to)  s  describes  the  rate  at  which  the  co-phasal 

surfaces  propagate  through  the  dispersive  medium  [2],  Since  the  phase  of  a  spatially 
coherent  optical  field  can  only  be  measured  indirectly  [52],  this  velocity  measure  does 
not  have  any  separate,  measurable  physical  meaning  in  spite  of  the  fact  that  it  plays 
a  central  role  in  the  mathematical  description  of  pulse  dispersion,  as  described  by 
the  Fourier-La  place  integral  representation  given  in  equation  (I),  In  particular, 
the  phase  velocity  of  a  pulse  is  superluminal  (Le,  vp(oO  >  c)  when  the  input  pulse 
carrier  frequency  is  above  the  uppermost  absorption  band  of  a  Lorentz  model 
dielectric,  as  illustrated  in  figure  1  for  a  double  resonance  Lorentz  model  dielectric 
when  toc  >  to 3. 

The  classical  group  velocity  vg(a>)  =  fk^r1  describes  the  rate  at  which  the 

envelope  of  a  group  of  waves  travels  through  the  dispersive  medium  [1-4]. 
As  described  by  Rayleigh  [2],  a  group  of  waves  is  defined  as  moving  beats  following 
each  other  in  a  regular  pattern  as,  for  example,  that  obtained  from  the  coherent 


Figure  1.  Angular  frequency  dispersjon  of  the  relative  phase  velocity  Vp/e  (dot-dashed 
curve),  relative  group  velocity  vg/c  (dashed  curve)  and  relative  energy  velocity  (solid  curve) 
in  a  double  resonance  Lorentz  model  dielectric  with  infrared  (wq  =  1 .74  x  10]4  rs^1, 
bo  =  L22  x  10ld  rs”\  =  4.96  x  10n  rs_t)  and  near  ultraviolet  (^  =  9. 145  x  IQ1*  r  s~'( 
bj  =  6 .72  x  10Jsrs~\  Si  =  1.434  x  10,5rs"j  resonance  lines. 
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corresponds  to  a  generalized  Brillouin  precursor  (describing  the  low-frequency 
response  of  the  medium  below  the  lowermost  absorption  band),  and  the  field 
component  /),  if  it  is  present,  corresponds  to  a  generalized  middle  precursor 
(describing  the  intermediate-frequency  response  of  the  medium  below  the  uppermost 
absorption  band).  A  genera)  condition  for  the  appearance  of  the  middle  precursor  is 
given  in  [42],  Each  field  component  contains  a  gaussian  amplitude  factor,  the 

peak  amplitude  point  of  each  propagating  at  the  classical  group  velocity  evaluated  at 
the  instantaneous  oscillation  frequency  of  that  field  component  at  that  space-time 
point  [54,  55],  As  the  pulse  evolves,  the  instantaneous  oscillation  frequency  at  each 
peak  amplitude  point  changes,  the  change  being  most  rapid  during  the  initial  pulse 
evolution  and  in  the  anomalous  dispersion  regions  where  the  material  dispersion  is 
greatest.  If  the  input  pulse  carrier  frequency  a>e  is  situated  in  a  region  of  anomalous 
dispersion  where  the  classical  group  velocity  is  either  superluminal  or  negative 
(or  both),  then,  as  the  pulse  begins  to  evolve  with  increasing  propagation  distance,  its 
instantaneous  oscillation  frequency  at  that  peak  amplitude  point  rapidly  shifts  into 
a  region  of  normal  dispersion  where  the  group  velocity  is  subluminal  and  approaches 
the  energy  transport  velocity  [55],  Hence,  any  superluminal  or  negative  movement  of 
the  peak  pulse  amplitude  point  occurs  in  the  initial  pulse  evolution  and  is  extremely 
short-lived.  Whether  or  not  this  motion  has  any  physical  significance  is  debatable  as 
peak  amplitude  points  in  the  dynamical  pulse  evolution  are  arguably  not  causally 
related  [56]  and  may  or  may  not  convey  any  transfer  of  information  [57-59], 
Nevertheless,  Sommerfeld’s  theorem  clearly  establishes  that  the  front  of  the  electro¬ 
magnetic  field  cannot  move  superluminally.  Since  the  field  at  some  propagation 
distance  can  always  be  separated  into  two  parts,  one  preceding  a  given  point  /0 
in  time  and  the  other  following  that  point  in  time,  then  for  all  larger  propaga¬ 
tion  distances,  Sommerfeld’s  theorem  states  that  the  field  part  following  that  initial 
point  in  time  cannot  move  ahead  of  the  propagated  space-time  point  0F  = 
(c/A:)(t  —  /o)  =  I  travelling  at  the  speed  of  light  c;  that  is,  peak  amplitude  points 
at  two  different  space-time  points  in  the  propagated  field  evolution  that  are  super¬ 
luminally  separated  are  not  causally  connected, 

It  is  clear  that  a  more  physically  meaningful  pulse  velocity  measure  needs  to  be 
considered  in  order  to  more  accurately  describe  the  complicated  pulse  evolution  that 
occurs  in  ultrashort  dispersive  pulse  dynamics.  One  possible  measure  is  given  by  the 
pulse  centrovelocity  [46] 


uCe  - 


lE2(r,t)  d  t 


£'(r.  /)dr 


(19) 


which  describes  the  evolution  of  the  temporal  centre  of  gravity  of  the  pulse  intensity. 
A  more  appropriate  velocity  measure  would  track  the  temporal  centroid  of  the 
Poynting  vector  of  the  pulse.  This  pulse  centroid  velocity  of  the  Poynting  vector  was 
first  introduced  by  Lisak  [47]  in  1976.  Recent  descriptions  [48-50]  of  its  properties 
have  established  its  efficacy  in  describing  the  evolution  of  the  pulse  velocity  with 
propagation  distance  in  a  Loren tz  model  dielectric, 

The  instantaneous  centroid  velocity  of  the  pulse  Poynting  vector  is  defined  as  [50] 

vci  =  lim  (Az/A (r))  (20) 
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Dynamical  Evolution  of  the  Brillouin  Precursor  in 
Rocard-Powles-Debye  Model  Dielectrics 

Kurt  Edmund  Oughstun 


Abstract — When  an  ultra  wide- band  electromagnetic  pulse  pen¬ 
etrates  into  a  causally  dispersive  dielectric,  the  interrelated  effects 
of  phase  dispersion  and  frequency  dependent  attenuation  alter 
the  pulse  In  a  fundamental  way  that  results  in  the  appearance 
of  so-called  precursor  fields.  For  a  Debye-type  dielectric,  the 
dynamical  field  evolution  is  dominated  by  the  Brillouin  precursor 
as  the  propagation  depth  typically  exceeds  a  single  penetration 
depth  at  the  carrier  frequency  of  the  input  pulse.  This  is  because 
the  peak  amplitude  in  the  Brillouin  precursor  decays  only  as  the 
square  root  of  the  inverse  of  the  propagation  distance.  This  nonex- 
ponential  decay  of  the  Brillouin  precursor  makes  it  ideally  suited 
for  remote  sensing.  Of  equal  importance  is  the  frequency  structure 
of  the  Brillouin  precursor.  Although  the  instantaneous  oscillation 
frequency  is  zero  at  the  peak  amplitude  point  oT  the  Brillouin 
precursor,  the  actual  oscillation  frequency  of  this  field  structure  is 
quite  different,  exhibiting  a  complicated  dependence  on  both  the 
material  dispersion  and  the  Input  pulse  characteristics.  Finally,  a 
Brillouin  pulse  is  defined  and  is  shown  to  possess  near  optimal  (if 
not  optimal)  penetration  into  a  given  Debye-type  dielectric. 

Index  Terms — Electromagnetic  propagation  in  absorbing  media, 
electromagnetic  propagation  in  dispersive  media,  electromagnetic 
transient  propagation,  ultra  wide  band  radar. 

I.  Introduction 

THE  dynamical  evolution  of  an  ultra  wide-band  electromag¬ 
netic  pulse  as  it  propagates  through  a  causally  dispersive 
dielectric  is  a  classical  problem  [  1  ]-[6]  in  electromagnetic  wave 
theory  with  considerable  current  importance  [7]-[!  1].  The  fre¬ 
quency  dependent  phase  and  attenuation  in  a  causal  medium 
are  interrelated  through  a  Hilbert  transform  [12].  Because  of 
this,  an  ultra  wide- band  pulse  undergoes  fundamental  changes 
as  it  propagates  through  a  dispersive  material.  Each  spectral 
component  travels  through  the  dispersive  medium  with  its  own 
phase  velocity  so  that  the  phasal  relationship  between  the  spec¬ 
tral  components  of  the  pulse  changes  with  propagation  distance, 
and  each  spectral  component  is  attenuated  at  its  own  rate  so 
that  the  relative  amplitudes  betwreen  the  spectral  components 
of  the  pulse  also  change  with  propagation  distance.  These  two 
interrelated  effects  result  in  a  complicated  dynamical  evolu¬ 
tion  of  the  propagated  field  [6]  that  is  accurately  described  by 
the  asymptotic  theory  [13  ]— f  17]  as  the  propagation  distance  ex¬ 
ceeds  a  value  set  by  the  material  absorption  depth  at  the  input 
pulse  carrier  frequency  [18].  For  an  ultrawide-band  pulse,  these 
combined  effects  manifest  themselves  through  the  formation 
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of  well-defined  precursor  fields  which  asymptotically  dominate 
the  dynamical  field  behavior  in  the  mature  dispersion  regime 
[19],  [20],  Notice  that  an  ultrawide- bund  pulse  need  not  be  ul- 
trashort;  for  example,  a  rectangular  envelope  modulated  pulse 
of  temporal  duration  T  >  0  has  a  speetmm  that  is  always  ul- 
trawide-band,  falling  off  as  \uf\~ 3  as  |w|  — >  oo,  while  the  pulse 
itself  may  not  be  ultrashort. 

The  precursor  fields  are  a  characteristic  of  the  material 
dispersion  [20],  the  input  pulse  merely  providing  the  requisite 
spectral  energy  in  the  appropriate  frequency  domain.  For  the 
Lorentz  model  [21]  of  resonance  polarization  phenomena,  used 
in  the  classical  theory  of  dispersive  pulse  propagation  [l]-[6], 
both  a  high-frequency  (above  resonance)  Sommerfeld  precursor 
and  a  low-frequency  (below  resonance)  Brillouin  precursor  are 
present  in  the  propagated  field  structure  when  the  input  pulse  is 
ultrawideband,  Additional  precursor  fields  may  also  exist  for  a 
multiple  resonance  Lorentz  model  dielectric  [22].  For  both  the 
Debye  model  [23]  of  orientational  polarization  phenomena  [24] 
and  the  Rocard-Powles  extension  [25]  of  the  Debye  model, 
only  the  Brillouin  precursor  field  is  present  in  the  propagated 
field  structure  [20].  Because  of  its  unique  nonexponentia!  peak 
decay,  the  Brillouin  precursor  has  direct  application  to  foliage 
and  ground  penetrating  radar,  remote  sensing  and  wireless 
communications  in  adverse  environments.  However,  its  effi¬ 
cacy  depends  upon  its  physical  frequency  structure  about  this 
peak  amplitude  point. 

For  carrier  frequencies  in  the  radio  spectrum  and  below  (fc  < 
300  GHz)  the  material  response  is  typically  dominated  by  ori¬ 
entational  polarization  effects  [24].  Above  this  frequency,  res¬ 
onance  polarization  effects  begin  to  dominate  the  material  dis¬ 
persion.  The  focus  of  this  paper  is  on  the  Brillouin  precursor  in 
Rocard-Powles-Debye  model  dielectrics  when  the  input  ultra-* 
wide-band  pulse  has  a  carrier  frequency  that  is  sufficiently  small 
that  resonance  polarization  effects  are  entirely  negligible. 

11,  Plane  Wave  Pulse  Propagation  in  Temporally 
Dispersive  Dielectrics 

The  angular  frequency  domain  form  of  Maxwell's  equations 
in  source- free  regions  of  a  temporally  dispersive  dielectric  is 
given  by  the  pair  of  equations  (with  t  = 

VxE(r  Tw)  = 

V  x  H(r,u>)  =  —  w)  (I) 

where  e(u?)  -  cr(u/)  +  tr,(w)  is  the  complex- valued  dielectric 
permittivity  whose  real  pan  cr(w)  =  3?{e(u>)}  and  imaginary 
part  r,(w)  =  3^e(u;)}  form  a  Hilbert  transform  pair  [12]  and 
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Fig.  1 .  Angular  frequency  dispersion  of  (he  real  {solid  curve)  and  imaginary 
(dashed  curve)  parts  of  the  complex  index  of  refraction  for  the  simple 
Rocard-Powles-Debye  model  of  triply-distilled  water 


presented  in  this  paper.  The  branch  points  of  n(w)  and  <f>(u>,0) 
include  the  singularities  wpi  —  -i/ry,  uP2  =  -i/r  and  zeroes 
w.i  =  -l(rp  -  c*)/(2t*J,  uK2  -  -1{t-p  +  c1)j( 2t£),  where 
rp  =  r  +  Tj,  =  rry,  and  c,  =  -  4r£(l  +  a/e*). 

For  a  Debye -type  dielectric,  the  saddle  point  equation  yields 
[20]  just  a  near  saddle  point  solution  in  the  low-frequency  do- 
main  about  the  origin.  For  |oj|  <  ju^*]  the  complex  index  of 
refraction  (10)  may  be  approximated  by  the  quadratic  expres¬ 
sion 


-Sr 


4f  -r* 


-(too  +  3£j)  -  1 


ft  flTfY 

U2  +  l-~h}  (11) 


2#o 


where  6q  =  n(0)  =  e lJ2 ,  With  this  substitution,  the  saddle 
point  equation  yields*  for  B  >  $0  -  k2/ the  approximate 
near  saddle  point  location1 


(12) 


with  k  =  arp/(2B0)  and  C  «  {aT^/(2Bo))[r2(€00  + 
3£J)/(4e-r^l)  —  1],  Numerical  results  show  that  this  saddle 
point  moves  down  the  imaginary  axis  as  $  increases  from  the 
value  Boo  =  crossing  the  origin  at  B  s=  00  and  then  ap¬ 
proaching  the  branch  point  singularity  wp2  —  -i/r  as  B  oo, 
as  described  by  (12), 

'This  approximate  solution  is  different  from  the  one  presented  in  [26,  eq. 
<21  >]  which  is  accurate  only  over  a  small  space-lime  interval  about  the  value 
H/tP  that  is  defined  in  that  paper 


For  B  >  Boo,  the  asymptotic  description  of  the  propagated 
field  in  a  Debye-type  dielectric  may  be  expressed  cither  in  the 
form  [20] 


A(z,t)  ^  i4fl(z,f)  +  Ae(M)  (13) 

as  A z  — *  oo,  or  else  in  a  form  that  is  a  superposition  of  expres¬ 
sions  of  the  form  given  in  (13).  For  example,  an  input  rectan¬ 
gular  envelope  pulse  of  temporal  duration  T  may  be  expressed 
as  the  difference  between  two  Heaviside  unit-step-function  sig¬ 
nals  displaced  in  time  by  T ,  each  propagated  signal  being  de¬ 
scribed  by  the  asymptotic  expression  given  in  (13),  the  first 
being  referred  to  as  the  leading-edge  signal  and  the  latter  as  the 
trailing -edge  signal.  Similar  representations  hold  for  other  ini¬ 
tially  symmetric  pulse  shapes.  The  field  quantity  Ab{z 1 1)  is  due 
to  the  asymptotic  contribution  from  the  near  saddle  point  and  is 
referred  to  as  the  Brillouin  precursor  [20],  The  field  contribution 
Ac(z ,  t)  is  due  to  the  pole  contributions  (if  any)  and  is  referred 
to  as  the  signal  component  [1  }-[*>]■ 

The  asymptotic  description  of  the  Brillouin  precursor  in  a 
Rocard-Powles-Debye  model  dielectric  is  obtained  through  a 
direct  application  [6]  of  Olver’s  theorem  [28]  to  the  contour 
integral  obtained  from  (8)  by  deforming  the  contour  C  to  an 
Giver-type  path  through  the  near  saddle  point,  with  the  result 
shown  in  (14)  at  the  bottom  of  the  page,  as  Az  — >  oo  with  B  > 
£oo  ,  where  w/v  (0)  denotes  the  near  saddle  point  location  and 
where  =  d2</>(w,B)/di J2  =  t[2n'(u;)  +  um"{ii;)].  Unlike 

that  for  a  Lorentz  model  dielectric  where  there  are  two  neigh¬ 
boring  near  saddle  points  that  coalesce  into  a  single  second- 
order  saddle  point,  thereby  requiring  uniform  [6],  [14]  and  tran¬ 
sitional  [29]  asymptotic  expansion  techniques,  the  asymptotic 
expression  given  in  (14)  is  uniformly  valid  for  all  finite  B  > 
da2  provided  that  any  pole  singularities  of  the  spectral  func¬ 
tion  u(dn{B)  —  utc)  are  sufficiently  well-removed  from  the  near 
saddle  point  location.  In  that  case,  the  pole  contribution  Ac{z,  i) 
is  given  by  a  direct  application  of  the  residue  theorem. 

For  example,  for  a  Heaviside  unit-step-function  signal,  u(^- 
ojc)  —  i/{wc  —  u>)  so  that  there  is  a  simple  pole  singularity  sit¬ 
uated  along  the  positive  real  frequency  axis  at  the  input  carrier 
frequency.  If  denotes  the  space-lime  value  when  the  Olvcr- 
type  path  P{9)  crosses  the  pole,  then  [6],  £13] — [  15J,  shown  in 
(35)  at  the  bottom  of  the  following  page,  as  Az  — *  oo  with 
B  >  The  uniform  asymptotic  theory  [6],  [14],  [30]  de¬ 
scribing  the  interaction  between  a  saddle  point  and  a  simple  pole 
singularity  of  the  integrand  identifies  0*  as  the  space- time  value 
B  at  which  the  steepest  descent  path  through  the  near  saddle 
point  crosses  the  simple  pole  singularity  at  —  u)c* 

IV.  Dynamical  Evolution  of  the  Brillouin  Precursor 

The  dynamical  structure  of  the  Brillouin  precursor  in  a  single 
relaxation  time  Rocard-Powles-Debye  model  dielectric  when 
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Fig.  3.  Dynamical  field  evolution  of  an  input  unit  amplitude  ten  cycle 
rectangular  envelope  pulse  with  fe  =  3  GHz  carrier  frequency  at  one.  three, 
five,  seven  and  nine  absorption  depths  in  the  simple  Rocard-Powles-Debye 
model  of  triply-di  stilted  water. 


Fig,  5.  Peak  amplitude  attenuation  as  a  function  of  the  relative  propagation 
distance  for  input  unit  amplitude  single  cycle  rectangular  envelope  pulses  with 
carrier  frequencies  }t  =  0.1  GHz  1*  symbols),  fc  =  1  GHz  (o  symbols),  and 
ft  =  10  GHz  {+  symbols).  The  solid  curve  describes  the  pure  exponential 
attenuation  given  by  exp{-A2/zd). 
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Fig,  4,  Propagated  pulse  sequence  of  a  unit  amplitude  rectangular  envelope 
single  cycle  pulse  with  /c  s=  1  GHz  earner  frequency  in  the  simple 
Rocard-Powles-Debye  model  of  triply-distilled  water 


when  the  propagation  distance  exceeds  a  single  absorption 
depth  at  the  pulse  carrier  frequency.  This  is  illustrated  in  the 
sequence  of  graphs  presented  in  Fig.  3  for  an  input  ten  cycle 
rectangular  envelope  pulse  with  1  GHz  carrier  frequency  and 
T  —  10  ns  initial  pulsewidth  at  one,  three,  five,  seven,  and 
nine  absorption  depths  in  water.  Notice  that  the  leading  and 
trai ling-edge  Brillouin  precursors  persist  long  after  the  1  GHz 
signal  has  been  significantly  attenuated  by  the  medium.  Al¬ 
though  these  two  Brillouin  precursors  penetrate  very  far  into 
the  material,  they  only  carry  a  small  fraction  of  the  initial  pulse 
energy  in  the  particular  case  illustrated  here.  The  leading -edge 
Brillouin  precursor  is  essentially  a  remnant  of  the  first  half  cycle 
of  the  initial  pulse  while  the  trai  ling-edge  Brillouin  precursor 
is  a  remnant  of  the  last  half  cycle.  The  input  pulse  energy 
available  to  the  leading  and  trail ing-edge  Brillouin  precursors 


is  then  limited  to  that  contained  in  a  single  cycle  of  the  input 
rectangular  envelope  pulse.  For  a  ten  cycle  pulse  as  illustrated 
in  Fig,  3*  this  means  that,  at  most,  only  ten  percent  of  the  input 
pulse  energy  is  available  to  this  precursor  pair  [33). 

A  more  efficient  way  to  generate  a  Brillouin  precursor  pair 
in  a  dispersive  material  is  with  a  single  cycle  pulse.  The  pulse 
sequence  presented  in  Fig.  4  illustrates  the  dynamical  pulse 
evolution  as  a  unit  amplitude,  rectangular  envelope  single 
cycle  pulse  with  fc  =  1  GHz  penetrates  into  water  with 
A z/zd  —  0, 1, 2, . . . ,  10.  The  evolution  of  the  pulse  into  a 
pair  of  leading  and  trad  ing-edge  Brillouin  precursors  is  clearly 
evident  as  the  propagation  distance  exceeds  a  single  absorp¬ 
tion  depth  (A zjzd  >  1)  and  the  peak  amplitude  attenuation 
goes  from  exponential  decay  to  the  (A^)*1^2  algebraic  decay 
described  in  (16).  This  transition  is  illustrated  in  Fig.  5  which 
presents  the  peak  amplitude  decay  as  a  function  of  the  relative 
propagation  distance  A z/zd.  The  solid  curve  in  the  figure 
describes  the  behavior  of  the  function  exp(-Az/zj)  which 
represents  pure  exponential  decay.  The  numerically  determined 
peak  amplitude  decay  for  three  different  input  single  cycle 
pulses  with  carrier  frequencies  /c  =  0,1  GHz,  fc  =  1  GHz, 
and  fc  =  10  GHz  is  presented  in  this  figure  by  the  ot  and 
4-  symbols,  respectively,  each  data  set  connected  by  a  cubic 
spline  fit 

The  temporal  width  of  the  leading-edge  Brillouin  precursor 
as  a  function  of  the  propagation  distance  is  illustrated  in  Fig.  6. 
Part  (a)  of  the  figure  is  plotted  in  terms  of  the  absolute  propaga¬ 
tion  distance  in  meters  while  (b)  is  plotted  in  terms  of  the  relative 
propagation  distance  A zjzd.  The  dependence  of  the  absorption 
depth  zd  —  a"1  (wc)  on  the  carrier  frequency  of  the  pulse  is  re¬ 
flected  in  the  individual  curves  appearing  in  Fig,  6<b),  The  solid 
curves  in  the  figure  describe  the  asymptotic  result  given  in  (20) 
and  the  three  sets  of  data  points  present  numerical  results  for  the 
fc  =  0,1  GHz,  /c  —  1  GHz,  and  fc  —  10  GHz  single  cycle 
pulse  cases,  each  data  set  connected  by  a  cubic  spline  fit.  Notice 
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Fig.  9.  Temporal  structure  of  the  Brill  cm  in  pulse  BPl  with  time  delay  T  s= 
l/{2/c)  for  /e  =  ]  GHz.  The  separate  leading  and  trailing-edge  precursor 
components  are  illustrated  in  part  (a)  and  their  superposition  is  given  in  pan  (b). 
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Fig.  10.  Propagated  pulse  sequence  for  the  Brillouin  pulse  BPl  with  delay  time 
T  =  1  /  (2/c)  for  fc  =  1  GHz  in  the  simple  Rocard- Pewits- Debye  model  of 
Iriply-distiiled  water. 


penetrate  a  finite  distance  into  a  given  dispersive  dielectric.  The 
results  presented  in  Figs.  4  and  5  indicate  that  the  pulse  that 
will  provide  near-optimal,  if  not  indeed  optimal,  penetration  is 
comprised  of  a  pair  of  Brillouin  precursor  structures  with  the 
second  precursor  delayed  in  time  and  ?r  phase  shifted  from  the 
first.  This  so-called  Brillouin  pulse  is  obtained  from  (14)  with 
A z  =  Zd  =  er_1^2(wc)  in  the  exponential,  the  other  factors  not 
appearing  in  the  exponential  set  equal  to  unity,  and  is  given  by 


fBp{t)  =  exp 


0(^jv(0M) 

.  WcTli(u>c)  _ 


—  exp 


uJcni(vc)  . 


(22) 


where  9r  =  8  —  cT/zd  with  T  >  0  describing  the  fixed 
time  delay  between  the  leading  and  trail  ing-edge  Brillouin  pre¬ 
cursors,  If  T  is  chosen  too  small  then  there  will  be  significant 
destructive  interference  between  the  leading  and  trail  ing-edges 
and  the  pulse  will  be  rapidly  extinguished.  For  practical  rea¬ 
sons,  2 T  should  be  chosen  near  to  the  inverse  of  the  operating 
frequency  fc  of  the  antenna  used  to  radiate  this  B  rillouin  pulse. 

With  T  —  I/(2/c)  the  input  Brillouin  pulse  is  approximately 
a  single  cycle  pulse  with  effective  oscillation  frequency  equal 
to  fc-  The  input  Brillouin  pulse  when  fc  =  1  GHz  is  depicted 
in  Fig,  9;  part  (a)  of  the  figure  shows  the  separate  leading  and 
trailing -edge  Brillouin  precursor  structures  while  part  (b)  of  the 
figure  shows  the  final  pulse  obtained  from  the  superposition  of 
these  two  parts,  as  described  by  (22),  The  initial  rise  and  fall 
time  for  this  pulse  is  ^  0.6  ns.  The  dynamical  evolution  of  this 
input  Brillouin  pulse  in  triply-distilled  water  is  illustrated  by  the 
pulse  sequence  given  in  Fig.  10  with  A z/zd  =  0, 1,2, ,  10. 
Comparison  with  the  pulse  sequence  depicted  in  Fig.  4  for  a 
1  GHz  single  cycle  rectangular  envelope  pulse  shows  that  the 
Brillouin  pulse  decays  much  slower  with  propagation  distance. 
Improved  results  are  obtained  when  the  delay  is  doubled  to  the 
value  T  =  l//c.  In  this  case  there  is  a  noticeable  “dead-time" 
between  the  leading  and  traiiing-edge  Brillouin  precursor  struc¬ 
tures  which  decreases  the  effects  of  destructive  interference  be¬ 
tween  these  two  components  of  the  Brillouin  pulse,  resulting 


in  improved  penetration  into  the  dispersive,  absorptive  mate¬ 
rial  However,  this  destructive  interference  can  never  be  com¬ 
pletely  eliminated  for  all  propagation  distances  since  the  time 
delay  between  the  peak  amplitude  points  for  the  leading  and 
traiiing-edge  Brillouin  precursors  decreases  with  the  inverse  of 
the  propagation  distance  [6],  [15]*  Nevertheless,  it  can  be  ef¬ 
fectively  eliminated  over  a  given  finite  propagation  distance  by 
choosing  the  time  delay  T  sufficiently  large.  The  tradeoff  in 
doing  this  is  to  decrease  the  effective  oscillation  frequency  of 
the  radiated  pulse. 

The  numerically  determined  peak  amplitude  decay  with  rel¬ 
ative  propagation  distance  A  z/zd  is  presented  in  Fig,  II.  The 
lower  solid  curve  depicts  the  exponential  attenuation  described 
by  the  function  exp{-Az/zd),  and  the  lower  dashed  curve  de¬ 
scribes  the  peak  amplitude  decay  for  a  single  cycle  rectangular 
envelope  pulse  with  fc  =  1  GHz.  Notice  that  the  departure  from 
pure  exponential  attenuation  occurs  when  A zjzd  ~  0.6  as  the 
leading  and  traiiing-edge  Brillouin  precursors  begin  to  emerge 
from  the  pulse.  The  dashed  curve  labeled  BP  I  describes  the  peak 
amplitude  decay  for  the  Brillouin  pulse  with  T  =  1/(2 fc)t 
BP2  describes  that  for  the  Brillouin  pulse  with  T  -  l//ct 
and  BPS  describes  that  for  T  =  3/(2/*),  There  is  no  notice¬ 
able  improvement  in  the  peak  amplitude  decay  as  the  delay 
time  T  is  increased  beyond  3/(2/r)  over  the  illustrated  range 
of  propagation  distances.  Notice  that  at  ten  absorption  depths, 
exp(-A zjzd)  -  exp(-10)  9*  4,54  x  10~\  the  peak  ampli¬ 
tude  of  the  single  cycle  pulse  is  0.07 1 8,  the  peak  amplitude  of 
the  Brillouin  pulse  BPl  is  0.2123,  the  peak  amplitude  of  the 
Brillouin  pulse  BP2  is  0.2943,  and  the  peak  amplitude  of  the 
Brillouin  pulse  BP3  is  0.3015,  over  three  orders  of  magnitude 
larger  than  that  expected  from  simple  exponential  attenuation. 

The  power  associated  with  the  observed  peak  amplitude 
decay  presented  in  Fig.  1 1  may  be  accurately  determined  [34] 
by  plotting  the  base  ten  logarithm  of  the  peak  amplitude  data 
versus  the  base  ten  logarithm  of  the  relative  propagation  dis¬ 
tance,  If  the  algebraic  relationship  between  these  two  quantities 
is  of  the  form  A p^k  =  B{Azfzd)p  where  B  is  a  constant,  then 
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seminal  analysis,  Brillouin  concluded  that11,12  "The  signal  velocity  does  not  differ  from  the  group 
velocity ,  except  in  the  region  of  anomalous  dispersion.  There  the  group  velocity  becomes  greater 
than  the  velocity  in  vacuum  if  the  reciprocal  c/U  <  I;  it  even  becomes  negative  ...Naturally,  the 
group  velocity  has  a  meaning  only  so  long  as  it  agrees  with  the  signal  velocity.  The  negative  parts 
of  the  group  velocity  have  no  physical  meaning. ..The  signal  velocity  is  always  less  than  or  at 
most  equal  to  the  velocity  of  light  in  vacuum”  This  research  then  established  the  asymptotic 
theory  of  pulse  propagation  in  dispersive,  absorptive  media.  An  essential  feature  of  this  approach 
is  its1  adherence  to  relativistic  causality  through  careful  treatment  of  the  dispersive  properties  of 
both  the  real  and  imaginary  parts  of  the  complex  index  of  refraction. 

At  approximately  the  same  time,  Havelock14 J 5  completed  his  research  on  wave  propagation  in 
dispersive  media  based  upon  Kelvin's  stationary  phase  method16.  It  appears  that  Havelock  was 
the  first  to  employ  the  Taylor  series  expansion  of  the  wave  number  (k  in  Havelock's  notation) 
about  a  given  wavenumber  value  Kq  that  the  spectrum  of  the  wave  group  is  clustered  about, 

referring  to  this  approach  as  the  group  method.  In  addition,  Havelock15  stated  that  "The  range  of 
integration  is  supposed  to  be  small  and  the  amplitude,  phase  and  velocity  of  the  members  of  the 
group  are  assumed  to  be  continuous ,  slowly  varying ,  functions  of  K.”  This  research  then 
established  the  group  velocity  method  for  dispersive  wave  propagation.  Since  the  method  of 
stationary  phase17  requires  that  the  wavenumber  be  real- valued,  this  method  cannot  properly  treat 
causally  dispersive,  attenuative  media.  Furthermore,  notice  that  Havelock’s  group  velocity 
method  is  a  significant  departure  from  Kelvin's  stationary  phase  method  with  regard  to  the 
wavenumber  value  Ko  about  which  the  Taylor  series  expansion  is  taken.  In  Kelvin's  method,  Ko  is 

the  stationary  phase  point  of  the  wavenumber  k  while  in  Havelock’s  method  describes  the 
wavenumber  value  about  which  the  wave  group  spectrum  is  peaked.  This  apparently  subtle 
change  in  the  value  of  Ko  results  in  significant  consequences  for  the  accuracy  of  the  resulting 
group  velocity  description. 

There  were  then  two  different  approaches  to  the  problem  of  dispersive  pulse  propagation:  the 
asymptotic  approach  (based  upon  Debye’s  method1"  of  steepest  descent)  which  provided  a  proper 
accounting  of  causality  but  was  considered  to  be  mathematically  unwieldy  without  any  simple, 
physical  interpretation,  and  Havelock’s  group  velocity  approximation  (based  upon  Havelock's 
reformulation  5  of  Kelvin’s  asymptotic  method15  of  stationary  phase)  which  violates  causality 
but  possesses  a  simple,  physically  appealing  interpretation.  It  is  interesting  to  note  that  both 
methods  are  based  upon  an  asymptotic  expansion  technique  but  with  two  very  different 
approaches,  the  method  of  stationary  phase  relying  upon  coherent  interference  and  the  method  of 
steepest  descent  relying  upon  attenuation. 

The  asymptotic  approach  was  revisited  by  Baerwald1*  in  1930  who  reconsidered  Brillouin ’s 
description1011  of  the  signal  velocity  in  causally  dispersive  systems,  and  also  by  Stratton14  in 
1941,  who  reformulated  the  problem  in  terms  of  the  Laplace  transform  and  derived  an  alternate 
contour  integral  representation  of  the  propagated  signal,  Stratton  appears  to  have  first  referred  to 
the  forerunners  described  by  Sommerfeld10  and  Brillouin11,12  as  precursors.  The  first  experimental 
measurement  of  the  signal  velocity  was  attempted  by  Shiren20  in  1962  using  pulsed  microwave 
ultrasonic  waves  within  a  narrow  absorption  band.  His  experimental  results20  were  * found  to  lie 
within  theoretical  limits  established  by  calculations  of  Brillouin  and  Baerwald *'  However,  a  more 
detailed  analysis  of  these  experimental  results  by  Weber  and  Trizna"3  indicated  that  the  velocity 
measured  by  Shiren  was  in  reality  that  for  the  first  precursor  and  not  the  signal.  Subsequent 
research  by  Handelsman  and  Bleistein22  in  1969  provided  a  uniform  asymptotic  description  of  the 
arrival  and  initial  evolution  of  the  signal  front.  The  first  experimental  measurements  of  the 
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Aco/(Or  « 1  is  satisfied.  The  frequency  dependence  of  the  wavenumber  may  then  be 
approximated  by  the  first  few  terms  of  its  Taylor  series  expansion  about  the  characteristic  pulse 
frequency  Q)e  with  the  assumption35,36,50  that  improved  accuracy  can  always  be  obtained  through 

the  inclusion  of  higher-order  terms;  this  assumption  has  been  proven  incorrect48’49,  optimal  results 
being  obtained  using  either  the  quadratic  or  the  cubic  dispersion  approximation  of  the 
wavenumber. 

Because  of  the  slowly-varying  envelope  approximation  together  with  the  neglect  of  the  frequency 
dispersion  of  the  material  attenuation,  the  group  velocity  approximation  is  invalid  in  the  ultrashort 
pulse  regime  in  a  causally  dispersive  material  or  system,  its’  accuracy  decreasing  as  the 
propagation  distance  Az  increases.  This  is  in  contrast  with  the  modem  asymptotic  description 
whose  accuracy  increases  in  the  sense  of  Poincare'7  as  the  propagation  distance  increases.  There 
is  then  a  critical  propagation  distance  ze  >  0  such  that  the  group  velocity  description  using  either 
the  quadratic  or  cubic  dispersion  approximation  provides  an  accurate  description  of  the  pulse 
dynamics  when  0  £  Az  <  zf,  the  accuracy  increasing  as  Az  — »0,  while  the  modem  asymptotic 

theory  provides  an  accurate  description  when  Az  >  zc,  the  accuracy  increasing  as  A z  — >  ».  This 
critical  distance  zc  depends  upon  both  the  dispersive  material  and  the  input  pulse  characteristics 
including  the  pulse  shape,  temporal  width  and  characteristic  angular  frequency  a)c .  For  example, 
z,  =  »  for  the  trivial  case  of  vacuum  for  all  pulse  shapes,  whereas  zr  -  zd  for  an  ultrashort, 
ultrawideband  pulse  in  a  causally  dispersive  dielectric  with  e~'  penetration  depth  zd  at  the 
characteristic  frequency  (Oe  of  the  input  pulse. 

In  spite  of  these  results,  the  group  velocity  approximation  remains  central  to  the  description  of 
ultrashort  pulse  dynamics  in  both  linear  and  nonlinear  optics  with  little  regard  to  its  domain  of 
validity.  This  is  seemingly  supported  by  the  apparent  agreement  between  experimental 
measurements  and  results  predicted  by  the  group  velocity  approximation.  Herein  lies  the  central 
controversy  considered  in  this  paper.  Related  to  this  is  the  controversy  regarding  the  possibility  of 
superluminal  pulse  velocities  since  the  group  velocity  can  assume  any  value  between  -»  and 
-h»  in  a  region  of  anomalous  dispersion. 

20.2.  Integral  Representation  of  the  Propagated  Pulse  and  Causality 

The  propagated  plane  wave,  pulsed  optical  field  A(z,t)  that  results  from  the  initial  pulse 
A{zo,/)  =  /(f)  at  the  plane  z  =  Zq  is  given  by  the  Fourier-Laplace  integral  representation 46 

=  (20.1) 

for  all  Az>0.  Here  f{(o)  is  the  temporal  angular  frequency  spectrum  of  the  initial  pulse 
function  /(r),  C  denotes  the  contour  of  integration  co  =  co'  +  ia  where  0)'  =  9?{a)}  ranges  from 
negative  to  positive  infinity,  and  a  is  a  constant  greater  than  the  abscissa  of  absolute  convergence 
for  /(f).  The  spectrum  A(z,(o)  of  the  optical  field  A(z,t)  satisfies  the  Helmholtz  equation 

[v2  +P(fo)]/4(z,o?)  =  0,  (20.2) 


where 
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20.3.  Havelock’s  Classical  Group  Velocity  Approximation 

The  group  velocity  approximation  is  a  hybrid  time  and  frequency  domain  representation50  in 
which  the  temporal  pulse  behavior  is  separated  into  the  product  of  a  slowly- varying  envelope 
function  and  an  exponential  phase  term  whose  angular  frequency  is  centered  about  some  fixed 
characteristic  frequency  of  the  initial  pulse.  Consider  then  the  specific  form  of  the  initial  pulse  at 
the  plane  Z  =  Zo  that  is  given  by  /(f)  =  u(r)sin(<yl.r  +  ^)  with  envelope  u(t)  and  constant  carrier 
frequency  OJr.  The  propagated  plane  wave  pulse  is  then  given  by  the  Fourier-Laplace  integral 
representation40 


A(z,t)  =  ^  Jr“(C0-a,r)e  W,'^£y| 


(20,7) 


for  all  Az  £0,  where  t (/  =  0,tt/2  for  either  a  cosine  or  sine  wave  carrier,  respectively.  Here 
ii(tu)  is  the  temporal  angular  frequency  spectrum  of  the  initial  pulse  envelope  function  u{t).  In 
the  slowly-varying  envelope  approximation,  the  envelope  function  «(/)  is  assumed  to  be  slowly- 
varying  on  the  time  scale  Aic  -  \j(Oc ,  which  is  equivalent51  to  the  quasimonochromatic 
approximation  that  the  spectral  bandwidth  A (O  of  u(o>)  is  sufficiently  narrow  that  the  inequality 
Aco/ojc  «1  is  satisfied.  The  complex  wavenumber  k(a>)  is  then  expanded  in  a  Taylor  series 


about  the  carrier  frequency  (Oc  with  the  assumption  '  that  this  series  may  be  truncated  after  a 
few  terms  with  some  undefined  error.  It  is  typically  assumed  that  the  attenuation  coefficient 
a{a))  =  3j£(a))J  is  sufficiently  small  that  it’s  frequency  dispersion  is  entirely  negligible  in 

comparison  to  that  for  the  propagation  factor  f3(a))  =  so  that  a{(o)  ~  a(a)r);  this  is 

entirely  compatible  with  the  stationary  phase  foundation14"10  of  the  group  velocity  description 
which  requires  that  the  wavenumber  be  real-valued  since  this  frequency  independent  attenuation 
factor  may  then  be  taken  outside  of  the  integration.  In  addition,  the  propagation  factor  is 

typically  represented  by  the  quadratic  dispersion  approximation 


p((Q) »  p((0f )+ ' \<oc )(oj  -  (oc )  +  —  p{2) (toc )(o>  -£U,.)\ 


(20.8) 


where  pt'j\(o)  =  dJ p(a))/do)J .  The  coefficient  /f(1,(cor)  is  the  inverse  of  the  group  velocity 
evaluated  at  the  carrier  frequency,  while  the  coefficient  )  describes  the  so-called  group 

velocity  dispersion so.  With  this  substitution,  Eq.  (20.7)  becomes 


A(z,r) ! 


f-a(<ur)Ai 


Si 


i[0(aic  )Ai-ov-r-3K/4] 


|  «(t')exP 


-I 


(fl,(o>r)Az  +  //-r)2 

2p(2](o)c)Az 


dt’  >,  (20.9) 


and  the  pulse  phase  propagates  through  the  dispersive  medium  at  like  phase  velocity 

CO 


vJu))m 


m' 


(20.10) 


while  the  pulse  envelope  propagates  through  the  dispersive  medium  at  the  group  velocity 
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field.  A  complete  understanding  of  these  saddle  point  dynamics  together  with  the  manner  in 
which  they  interact  with  the  initial  pulse  spectrum  provides  a  detailed,  accurate  description  of  the 
entire  dynamical  evolution  of  the  propagated  pulse  in  the  dispersive,  absorptive  medium  for  all 
A z>z(,  the  accuracy  of  this  approximation  increasing  in  the  sense  of  Poincare1'  as  the 
propagation  distance  increases  above  some  critical  propagation  distance  ze  >0. 

The  set  of  saddle  points  of  the  complex  phase  function  <p((0,B)  =  i£i)[n(£0)  —  0]  is  determined  by 
the  condition  that  0(tt),0)  be  stationary  at  a  saddle  point,  in  which  case  $'((0,6)  =  0,  where  the 
prime  denotes  differentiation  with  respect  to  £0,  so  that 

n(Q))  +  0)n'(0})  =  d  .  (20.14) 

The  solutions  of  this  saddle  point  equation  then  give  the  desired  saddle  point  locations  in  the 
complex  di-plane  as  a  function  of  the  space-time  parameter  0  =  ct/Az.  The  saddle  points  will 
then  evolve  with  time  at  any  fixed  propagation  distance  Az .  Because  of  the  general  symmetry 
relations43,4*  «(—££>)  -  n'(oi")  and  </>(-££),  0)  —  that  are  satisfied  by  a  causal  medium,  if 

£0,(0)  is  saddle  point  of  <f>( £0,0),  then  so  also  is  -£O*(0). 


An  appreciation  of  the  physical  significance  of  the  saddle  points  can  be  obtained  from  the  relation 
(Az/c)<f>((0,9)  =  tm£ti)Az  —  £Uf)  for  the  complex  phase  function.  Upon  differentiating  this 


expression  with  respect  to  (0,  one  obtains 


(A z/c)<p'((0,6)  =  j’^8A(£0)/8q)JAz  —  /J ,  Since 


(f>'((0,6)  =  0  at  each  saddle  point  £0,  (0)  of  <f> ,  then 


A z 

t 


1 

(3fc(£O)/0£o) 

\  /  at  =£0 


(20.15) 


and  the  complex  group  velocity  is  real-valued  at  the  saddle  points. 

With  the  saddle  point  locations  known  for  0  >  1 ,  the  asymptotic  analysis  then  proceeds  by 
expressing  the  integral  representation  given  in  Eq,  (20.5)  in  terms  of  an  integral  7(z,0)  with  the 
same  integrand  but  with  a  new  contour  of  integration  P(6)  to  which  the  original  contour  C  may 
be  deformed43"46.  By  Cauchy’s  residue  theorem,  the  integral  representation  (20.5)  of  A(z,t)  and 
the  contour  integral  7(z,0)  are  related  by 


where 


A(z,r)  =  7(z,0)  —  9?{2ffl‘A(0)} 


(20.16) 


A(0)  =  X^|^/(£u)exp[(Az/c)0(£O,0)]| 


(20.17) 


is  the  sum  of  the  residues  of  the  poles  that  were  crossed  in  the  deformation  from  C  to  ^(0),  and 
where 
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values  of  z>  Zq  for  sufficiently  large  t ,  As  a  result,  this  contribution  to  the  asymptotic  behavior  of 

the  propagated  field  describes  the  steady-state  behavior  of  the  signal.  The  arrival  of  this  signal 
contribution  is  determined  by  the  dynamics  of  that  dominant  saddle  point  that  becomes 
exponentially  negligible  in  comparison  to  the  pole  contribution.  A  detailed  knowledge  of  the 
saddle  point  dynamics  for  a  given  dispersive  material  is  then  seen  to  be  a  critical  ingredient  for  a 
detailed  description  of  dispersive  pulse  propagation  in  that  material,  not  just  for  the  transient  field 
behavior  described  by  Eq.  (20,19),  but  also  for  the  steady-state  behavior  described  in  Eq.  (20,20). 


In  a  multiple  resonance  Lorentz  model56  dielectric  the  complex  index  of  refraction  is  described 
by46 


»(")  = 


N 

1-1 


*>1 


j- o  ^  ^2  j  +  2/<52  jO)  j 


(20,21) 


where  C02J  is  the  undamped  resonance  frequency,  b2j  is  the  plasma  frequency,  and  S2J  the 


phenomenological  damping  constant  for  the  (2j/h  resonance  line  of  the  dielectric  material.  This 
causal  model52  provides  an  accurate  description  of  both  the  normal  and  anomalous  dispersion 
phenomena  observed  in  homogeneous,  isotropic,  locally  linear  optical  materials.  The  regions  of 

anomalous  dispersion  approximately  extend  over  each  frequency  domain  ^GJ2j,£02;+1  J,  where 
2j+i  =  +  -  In  this  case  the  saddle  point  equation  (20,14)  has  at  least  two  sets  of  saddle 


CO 


points  that  are  symmetrically  situated  about  the  imaginary  axis.  One  pair  of  saddle  points  (the 
distant  saddle  points  n'12’43,46)  evolve  in  the  high-frequency  region  \o>\  >  ti)2v+i  the  complex  Q)~ 
plane  above  the  uppermost  absorption  band  of  the  material,  while  another  pair  of  saddle  points 
(the  near  saddle  poults11*12,43’46)  evolve  in  the  low-frequency  region  |o>|  <  CO0  of  the  complex  (O- 
plane  below  the  lowermost  absorption  band  of  the  dielectric.  If  the  dielectric  material  is  described 
by  multiple  resonance  lines,  then  additional  middle  saddle  points  will  appear  in  the  region 
M  <  (01N  below  the  uppermost  absorption  of  the  dielectric.  The  asymptotic  description  of  the 
propagated  pulse  may  then  be  expressed  either  in  the  form 


A{z,t)  ~  As(z,t)  +  Am{z,t)  +  AB(z,t)  +  Ac(zj)  (20.22) 


as  A z  — >  or  by  an  expression  that  is  a  superposition  of  expressions  of  the  form  given  in  Eq, 

(20,22). 

Here  As(z,f)  denotes  the  contribution  from  the  distant  saddle  points  with  nonuniform  asymptotic 
approximation  given  by  Eq.  (20.19)  for  6>  1  and  is  referred  to  as  the  first  or  Sommerfeld 
precursor.  This  non  uni  form  approximation  breaks  down  at  8  —  1  when  the  distant  saddle  points 
are  at  infinity.  The  uniform  asymptotic  description2245,46  of  the  Sommerfeld  precursor,  uniformly 
valid  in  the  space-time  parameter  6  —  ctjAz  for  all  1,  must  then  be  used  in  place  of  Eq, 
(20.19)  for  the  initial  pulse  evolution.  The  instantaneous  angular  oscillation  frequency  of  the 
Sommerfeld  precursor  is  approximately  given  by  the  real  part  of  the  distant  saddle  point 
location43,46  in  the  right-half  of  the  complex  a; -plane.  The  Sommerfeld  precursor  then  describes 
the  signal  front  which  arrives  at  8  =  1  with  infinite  oscillation  frequency  (but  zero  amplitude  for 
a  finite  energy  input  pulse)  and  consequently  propagates  at  the  speed  of  light  c  in  vacuum.  As  8 
increases  away  from  unity,  the  amplitude  of  the  Sommerfeld  precursor  rapidly  increases  to  a 
maximum  value  and  then  decreases  monotonically  for  all  larger  8  while  the  instantaneous 
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20.5,  Accuracy  of  the  Group  Velocity  Description  of  Ultrashort  Pulse  Dynamics 

The  accuracy  of  the  group  velocity  approximation  of  ultrashort  pulse  dispersion  is  now 
considered  in  order  to  establish  the  space- time  domain  over  which  this  approximate  description  is 
valid.  A  double  resonance  Lorentz  model  of  a  fluoride-type  glass  with  infrared 
(Ct)0  =  1.74  x  10i4  r/s,  b0  =  1.22  x  1014  r/s,  SQ  =  4.96x1  Onr/s)  and  near  ultraviolet 

(co2  =  9.145 x  1015 r/s,  b2  =  6.72xl015r/j,  52  =  1.434 x  10u r/s)  resonance  lines  is 

considered  with  complex  index  of  refraction  given  by  Eq.  (20.21)  with  N  —  2.  The  angular 
frequency  dispersion  of  the  real  and  imaginary  parts  of  the  complex  wavenumber 
£(m)  =  (o}/c)n((o)  for  this  double  resonance  Lorentz  model  dielectric  is  illustrated  in  Figure  1. 
The  upper  and  lower  solid  curves  in  each  part  of  the  figure  describe  the  exact  frequency 
dependence  of  P(o))  =  and  a(at)s sjfcfw)),  respectively,  while  the  dashed  red  curves 

describe  the  cubic  dispersion  approximation  when  (a)  o)c  =  a)mi0  =  1.615x  10i5 r/s  and  (b) 
(Oc  =  0.87o>2  =  8.0 x  1015 rjs.  The  cubic  dispersion  approximation  is  seen  to  provide  a  reasonably 
accurate  estimate  of  the  local  frequency  dispersion  of  the  propagation  factor  /3(co)  about  the 
carrier  frequency  within  the  passband  where  the  dispersion  is  normal  when  (Oe  =  £0^.  but  the 
accuracy  of  this  approximation  is  seen  to  decrease48,49  as  <oc  is  shifted  toward  either  absorption 
band  where  the  dispersion  becomes  anomalous.  The  inclusion  of  higher-order  terms  in  the  Taylor 
series  approximation  of  the  complex  wavenumber  only  serves  to  further  decrease  its  accuracy  in  a 
global  sense48,49.  The  cubic  dispersion  approximation  of  the  attenuation  coefficient  is  not  as 
accurate  as  that  for  the  propagation  factor  making  the  necessity  of  the  approximation 
a(a))  =  a((Oc)  used  in  the  group  velocity  description  all  the  more  important. 


Figure  1.  Angular  frequency  dependence  of  the  real  (upper  solid  curves)  and  imaginary  (lower 
solid  curves)  parts  of  the  complex  wavenumber  for  a  double  resonance  Lorentz  model  of  a 
fluoride-type  glass  with  infrared  and  near  ultraviolet  resonance  lines.  In  pan  (a)  the  dashed 
curves  describe  the  cubic  dispersion  approximation  about  the  minimum  dispersion  point  in  the 
passband  between  the  absorption  bands  and  in  part  (b)  they  describe  that  approximation  about 
the  upper  resonance  frequency. 

Because  of  its  central  importance  in  ultrashort  optical  pulse  technology,  a  unit  amplitude  gaussian 
envelope  pulse  f(t)  =  w(r)sin(ft)cf  +  nf  2)  is  considered,  where 
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superluminal  and  negative.  Part  (a)  of  the  figure  illustrates  the  approximate  group  velocity  and 
actual  pulses  at  one  absorption  depth  A z-zd  in  the  medium,  where  zd  =  a~[(o)c)\  the  actual 
pulse  is  clearly  not  travelling  at  a  rate  given  by  the  classical  group  velocity  v^{cur)  =  l//J'(£Ur)( 

Part  (b)  of  the  figure  illustrates  the  actual  and  approximate  pulses  when  the  approximate  group 
velocity  pulse  (dashed  curve)  has  been  temporally  shifted  such  that  its  peak  amplitude  point  is 
coincident  with  that  for  the  actual  pulse  (solid  curve).  The  error  between  these  two  pulses  is  then 
computed  in  two  different  ways.  The  first  error  measure  (errorl)  is  given  by  the  integral  of  the 
square  of  the  difference  between  the  two  aligned  pulses.  This  error  then  measures  both  shape  and 
energy  differences.  The  second  error  measure  (error2)  is  obtained  by  first  renormalizing  both 
pulses  by  the  square  root  of  their  respective  pulse  energies  and  then  taking  the  integral  of  the 
square  of  their  differences.  This  error  then  measures  the  shape  difference  between  the  two  aligned 
pulses  at  a  given  fixed  propagation  distance. 


Figure  3.  Error  resulting  from  the  group  velocity  description  of  the  propagated  gaussian 
envelope  pulse  with  a  cubic  dispersion  approximation  of  the  propagation  factor  /?{&>)  in  a 
double  resonance  Lorentz  model  dielectric  as  function  of  the  relative  penetration  depth  Azfz^ 
when  the  input  pulse  carrier  frequency  a or  is  equal  to  the  angular  frequency  a)mm  at  the 

minimum  dispersion  point  in  the  passband  between  the  two  absorption  bands  for  different 
values  of  the  input  pulse  width  2  T . 

The  numerical  results  for  these  two  error  measures  are  depicted  in  Figure  3  as  a  function  of  the 
relative  propagation  distance  A zjzd  when  the  input  pulse  carrier  angular  frequency  COc  is  set 

equal  to  the  angular  frequency  0)^  at  the  minimum  dispersion  point  in  the  passband  between  the 
two  absorption  bands  of  the  double  resonance  Lorent2  model  dielectric.  The  results  are  presented 
for  four  different  values  of  the  input  gaussian  envelope  pulse  width  27\  The  results  clearly  show 
that,  as  expected,  the  error  decreases  with  increasing  initial  pulse  width  at  any  fixed  value  of  the 
propagation  distance.  The  first  error  measure  (enorl)  is  seen  in  part  (a)  of  the  figure  to  initially 
increase  with  increasing  propagation  distance,  reaching  a  peak  value  at  approximately  one 
absorption  depth  (A z(zd  1)*  and  then  decreasing  to  zero  as  the  relative  propagation  distance 
A zjzd  increases  above  unity.  This  behavior  is  due  to  the  fact  that  the  group  velocity  approximate 
pulse  amplitude  decays  exponentially  with  penetration  distance  at  a  faster  rate  than  does  the 
actual  ultrashort  pulse  so  that  the  observed  mono  tonic  decrease  in  error  with  increasing 
propagation  distance  A z/zd  >  1  primarily  describes  the  slow,  nonexponentiat  amplitude  decay 
with  propagation  distance  of  the  actual  pulse  in  the  lossy  dielectric.  This  difference  between  the 
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20.6.  The  Question  of  Superluminal  Pulse  Velocities 


A  number  of  velocity  measures  have  been  introduced  for  the  purpose  of  describing  the  rate  at 
which  some  particular  feature  of  a  pulse  travels  through  a  dispersive  material.  The  most  important 
of  these  are  the  phase3,  group1-4,30"32,  energy38,41,  signal5'12'18,4  ,  and  centroid61"65  velocities.  Each 

of  these  velocity  measures  in  free  space  equals  the  velocity  of  light  c  in  vacuum,  but  they  are  in 
general  different  in  a  causally  dispersive  material  such  as  that  described  by  the  Lorentz  model. 


The  phase  velocity  vp(£o)  =  (oj /3(cu)  describes  the  rate  at  which  the  co-phasal  surfaces  propagate 

through  the  dispersive  medium3  [cf.  Eqs.  (20.9)-(20.10)J.  Since  the  phase  of  a  spatially  coherent 
optical  field  can  only  be  measured  indirectly66,  this  velocity  measure  does  not  have  any  separate, 
measurable  physical  meaning  in  spite  of  the  fact  that  it  plays  a  central  role  in  the  mathematical 
description  of  pulse  dispersion,  as  described  by  Eq.  (20.1).  In  particular,  the  phase  velocity  of  a 
pulse  is  superluminal  [i.e.  v/,(tur)>c]  when  the  input  pulse  carrier  frequency  Q)e  is  above  the 

uppermost  absorption  band  of  a  Lorentz  model  dielectric,  as  illustrated  in  Fig.  5(a)  for  a  single 
resonance  Lorentz  model  dielectric  when  (Oe  >  €Ol .  For  an  ultrashort  pulse  with  above  resonance 
carrier  frequency  whose  temporal  energy  centroid  is  moving  subluminally,  the  phase  velocity  is 
then  seen  to  describe  the  motion  of  a  space-time  point  where  there  is  negligible  pulse  energy. 


Figure  5.  Frequency  dependence  of  (a)  the  relative  phase  velocity  (dashed  curve), 
group  velocity  (dash-dot  curve)  and  energy  velocity  (solid  curve)  and  (b)  the  relative 
energy  velocity  (solid  curve)  and  signal  velocity  (open  circles)  in  a  single  resonance 
Lorentz  model  dielectric. 


The  group  velocity  vJ(&>)  =  (3 ^{(6)fd(o)  '  describes  the  rate  at  which  the  envelope  of  a  group 

of  waves  travels  through  the  dispersive  medium1-4  [cf.  Eqs.  (20.9)  and  (20.11)].  As  described  by 
Rayleigh3,  a  group  of  waves  is  defined  as  moving  beats  following  each  other  in  a  regular  pattern 
as,  for  example,  that  obtained  from  the  coherent  superposition  of  two  monochromatic  waves  with 
slightly  different  amplitudes  and  frequencies.  Although  the  group  velocity  does  indeed  describe 
the  beat  velocity  of  such  an  infinite  wave  group,  its  extension  to  the  description  of  the  velocity  of 
an  ultrashort  pulse  in  a  causally  dispersive  medium  is  invalid.  This  is  readily  evident  in  Fig.  2(a) 
where  the  group  velocity  approximate  pulse  is  seen  to  be  moving  at  a  faster  rate  than  is  the  actual 
pulse.  This  example  illustrates  the  extreme  dispersion  case  when  0)c=Ci)2.  The  group  velocity 

value  v^cu,.)  is  then  seen  to  be  a  poor  measure  of  the  actual  pulse  velocity  when  the  material  loss 
is  not  negligible,  although  it  can  describe  the  initial  pulse  evolution  when  the  material  loss  is  near 
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For  the  numerical  gaussian  pulse  example  illustrated  in  Fig.  2,  the  peak  amplitude  point  of  the 
actual  pulse  is  moving  with  the  average  velocity  vfJJt  =  0.91c  and  the  peak  amplitude  point  of 

the  group  velocity  approximate  pulse  is  moving  with  the  average  velocity  =  —0.30c,  while 

the  phase  velocity  is  given  by  vp((Oc)  =  0.82c,  the  group  velocity  is  given  by  =  -0.41c, 

and  the  energy  velocity  is  given  by  v£(£Of)  =  0.18c  at  the  input  pulse  carrier  angular  frequency 

0)e  -  0t)2  =  9. 145  X  1015  rjs ,  The  peak  amplitude  values  reported  here  are  “time-of-flight”  values 
that  result  from  numerical  measurements  of  the  initial  and  final  pulse  positions.  This  average 
peak  amplitude  velocity  measure  changes  as  the  propagation  distance  into  the  dispersive 
dielectric  increases,  as  does  the  instantaneous  peak  amplitude  velocity.  The  instantaneous  peak 
amplitude  velocity  of  an  ultrashort  gaussian  pulse  has  been  shown59,60  to  evolve  with  increasing 
propagation  distance  along  the  group  velocity  curve  toward  the  energy  velocity  curve  as  the 
instantaneous  oscillation  frequency  at  the  peak  amplitude  point  shifts  away  from  the  region  of 
anomalous  dispersion  and  into  the  norma!  dispersion  region  either  above  or  below  that  absorption 
band.  It  is  then  not  surprising  that  the  numerical  velocity  values  given  above  are  all  significantly 
different  as  each  describes  a  different  feature  of  the  pulse  that  may  only  be  valid  either  in  the  limit 
of  vanishingly  small  propagation  distance  (the  group  velocity)  or  else  in  the  large  propagation 
distance  asymptotic  limit  (the  energy  velocity). 


Figure  6.  Frequency  dependence  of  the  relative  phase  velocity  (dashed  curves),  group  velocity 
(dash-dot  curves)  and  energy  velocity  (solid  curves)  for  (a)  the  double  resonace  Loren tz  model 
dielectric  considered  in  Fig.  I  and  (b)  the  same  double  resonance  Lorentz  model  dielectric 
when  each  of  the  phenomenological  damping  constants  have  been  reduced  by  a  factor  of  ten. 

It  is  clear  that  a  more  physically  meaningful  pulse  velocity  measure  needs  to  be  considered  in 
order  to  accurately  describe  the  complicated  pulse  evolution  that  occurs  in  ultrashort  dispersive 
pulse  dynamics.  One  possible  measure  is  given  by  the  pulse  centro velocity61 


v« s  m  j"  <irv,o*/£  ev.o*) 


(20.26) 


which  describes  the  temporal  center  of  gravity  of  the  pulse  intensity.  A  more  appropriate  velocity 
measure  would  track  the  temporal  centroid  of  the  Poynting  vector  of  the  pulse.  This  pulse 
centroid  velocity  of  the  Poynting  vector  was  first  introduced  by  Lisak6"  in  1976.  Recent 
descriptions63-65  of  its  properties  have  established  its  efficacy  in  describing  the  evolution  of  the 
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distances  to  its  asymptotic  behavior  for  large  propagation  distances.  The  minimum  point  in  this 
dip  occurs  at  a  propagation  distance  whose  value  increases  with  increasing  initial  pulse  width 
while  the  minimum  value  decreases  with  increasing  initial  pulse  width.  At  the  leading  edge  of 
each  dip  when  the  centroveiocity  rapidly  decreases,  the  pulse  begins  to  separate  into  a  pair  of 
middle  and  Brillouin  precursors  and  the  temporal  pulse  centroid  is  found  to  occur  at  a  space-time 
point  between  these  two  pulse  components  where  the  pulse  energy  is  minimal.  Finally,  notice  that 
the  centroveiocity  for  a  gaussian  pulse  is  subluminal  and  nonnegative  (i.e.,  0  <  vc/  <  c)  for  all 
propagation  distances  at  this  input  carrier  frequency. 

In  Fig.  7(b)  the  input  pulse  carrier  frequency  is  set  equal  to  the  upper  resonance  frequency  0)2  of 
the  double  resonance  Lorentz  model  dielectric  where  the  dispersion  is  anomalous.  The  classical 
group  velocity  limit  vg(eo()  =  -0.4076c  is  again  approached  at  a  sufficiently  small  propagation 

distance  in  the  limit  as  the  initial  pulse  width  increases  and  the  initial  pulse  spectrum  narrows 
about  the  carrier  frequency  C0c,  but  at  a  much  slower  rate  than  that  obtained  at  the  minimum 
dispersion  point.  However,  notice  that  if  the  initial  pulse  width  is  sufficiently  small  (as  it  is  for  the 

I37fs  pulse  case),  then  the  initial  pulse  spectrum  is  extremely  ultrawideband  such  that  the 
classical  group  velocity  limit  is  not  obtained  as  Az  — »  0  and  the  centroveiocity  remains  positive 
for  all  propagation  distances  Az  >  0 ,  In  the  opposite  limit  as  Az/zd  — >  00  the  centroid  velocity  is 

again  found  to  approach  the  velocity  vg  =  c/n(0)  that  the  peak  amplitude  point  in  the  Brillouin 
precursor  travels  at  through  the  dispersive  material.  For  a  sufficiently  long  initial  pulse  width,  the 
transition  of  the  ultrawideband  pulse  centroid  velocity  between  these  two  limits  is  marked  by  a 
rapid  decrease  in  centroveiocity  to  and  then  from  -H»  to  subluminal  values  before 
approaching  the  asymptotic  limit  vB  =  c/n(0)  set  by  the  peak  amplitude  point  of  the  Brillouin 
precursor.  The  discontinuous  jump  in  the  centroveiocity  from  -«  to  -H»  is  found  to  occur  at  a 
relative  propagation  distance  whose  value  increases  with  increasing  pulse  width  provided  that  the 
initial  pulse  spectrum  is  ultrawideband. 

In  spite  of  the  fact  that  the  instantaneous  centroveiocity  of  the  pulse  Poynting  vector  can  take  on 
both  negative  and  superluminal  values  for  sufficiently  small  relative  propagation  distances,  the 
pulse  itself  is  found  to  only  undergo  a  slight  change  in  shape.  There  isn’t  any  superluminal 
movement  of  the  pulse  when  the  instantaneous  pulse  centroveiocity  is  superluminal,  nor  is  there 
any  retrogression  in  position  when  the  pulse  centroveiocity  is  negative.  Sommerfeld's  theorem 
firmly  establishes  that  electromagnetic  field  energy  cannot  move  forward  of  any  space-time  point 
in  the  pulse  at  a  superluminal  rate.  Finally,  notice  that  this  seemingly  nonphysical  behavior  only 
occurs  in  the  immature  dispersion  regime  (0  ^  Az  <  zc)  where  the  group  velocity  approximation 

applies  in  the  limit  as  Az  — >  0 . 

20.7.  Conclusions 

The  analysis  and  numerical  results  presented  in  this  paper  have  established  the  following  results: 
(1)  The  group  velocity  approximation  is  valid  only  in  the  immature  dispersion  regime  0  <  Az  <  zr, 
its  accuracy  increasing  as  Az— >0.  The  asymptotic  description  is  valid  in  the  mature  dispersion 
regime  Az>zf,  its  accuracy  increasing  in  the  sense  of  Poincare1 7  as  Az  — >«.  The  critical 
distance  Zf  depends  upon  the  input  pulse  type  and  initial  pulse  length,  as  well  as  upon  the  input 
pulse  carrier  frequency  for  a  given  dispersive  material.  (2)  The  instantaneous  centroid  velocity  of 
the  pulse  Poynting  vector  is  a  convenient,  albeit  sometimes  misleading,  measure  of  the  pulse 
evolution  in  a  dispersive  medium  for  input  gaussian  envelope  pulses.  Although  this  velocity 
measure  can  take  on  both  negative  and  superluminal  values  for  relative  propagation  distances  in 
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The  classical  Lorentz  mode)  of  dielectric  dispersion  is  based  on  the  microscopic  Lorentz  force  relation  and  New¬ 
ton's  second  law  of  motion  for  an  ensemble  of  harmonically  bound  electrons.  The  magnetic  held  contribution  in 
the  Lorentz  force  relation  is  neglected  because  it  is  typically  small  in  comparison  with  the  electric  field  con¬ 
tribution.  Inclusion  of  this  term  leads  to  a  microscopic  polarization  density  that  contains  both  perpendicular 
and  parallel  components  relative  to  the  plane  wave  propagation  vector.  The  modified  parallel  and  perpendicu¬ 
lar  polarizabilities  are  both  nonlinear  in  the  local  electric  field  strength.  ©  2006  Optica)  Society  of  America 
OCIS  codes:  260.2030,  190.4400,  160.4760,  160.4330, 


1.  INTRODUCTION 

The  Lorentz  model of  resonant  dispersion  phenomena 
in  dielectric  materials  is  a  classical  model  of  central  im¬ 
portance  in  optics4,6  as  well  as  in  the  broader  discipline  of 
electromagnetics,6  It  is  a  causal  model6,7  that  describes 
both  normal  and  anomalous  dispersion  phenomena  from 
the  infrared  through  the  optical  regions  of  the  electro¬ 
magnetic  spectrum.  This  model  is  based  on  the  micro¬ 
scopic  Lorentz  force  relation 

1 

ffrtf)spOM)e(r,f)  +  — j(r»f)  x  h(rtt)t  (D 

where  p(r,f)  is  the  microscopic  change  density  and  j(r,£) 
=  p(r,f)dr/d/  the  microscopic  current  density  of  the 
charged  particle  at  the  space-time  point  (r,£)  and  where 
ftr,/)  is  the  microscopic  force  density  exerted  on  the 
charged  particle  by  the  electromagnetic  field  with  micro¬ 
scopic  electric  intensity  vector  e(r ,£)  and  magnetic  induc¬ 
tion  vector  b(r,0‘  The  quantity  *  appearing  in  the  double 
brackets  ||*||  in  any  equation  here  is  used  as  a  conversion 
factor  between  cgs  and  MKS  units,6'8  If  that  factor  is  in¬ 
cluded  in  that  equation,  then  the  equation  is  in  cgs  units, 
while  if  the  factor  is  omitted  from  that  equation,  then  it  is 
in  MKS  units.  If  no  such  factor  appears  in  a  given  equa¬ 
tion,  then  that  equation  is  correct  in  either  system  of 
units. 

The  magnetic  field  contribution  appearing  on  the  right- 
hand  side  of  Eq.  (1)  is  typically  assumed  to  be  negligible 
in  comparison  with  the  electric  field  contribution  for  suf¬ 
ficiently  small  field  strengths  and  values  of  the  relative 
velocity  vie  of  the  charged  particle,  where  u  =  |dr/dfj.  This 
assumption  then  leads  to  the  classicaJ  Lorentz  theory  of 
dielectric  dispersion  in  which  electric  field  effects  alone 
are  considered.  The  question  then  remains  as  to  what  ef¬ 
fects  are  introduced  by  the  inclusion  of  the  magnetic  field 
contribution  in  the  Lorentz  theory.  The  solution  of  this 
important  problem  has  been  partially  addressed9"11  for 


the  special  case  of  a  chiral  molecule,  where  it  has  been 
shown11  that  the  magnitude  of  the  magnetic  response  of 
the  induced  second-order  optical  activity  can  be  compa¬ 
rable  to  the  magnitude  of  the  electric  response  alone. 
However,  as  we  have  been  unable  to  find  any  previously 
published  general  solution  of  this  problem,  it  is  addressed 
in  this  paper. 


2.  MODIFIED  LORENTZ  MODEL  OF 
DIELECTRIC  DISPERSION 

With  the  complete  Lorentz  force  relation  given  in  Eq.  (I) 
as  the  driving  force,  the  equation  of  motion  of  a  harmoni¬ 
cally  bound  electron  is  given  by 


<u2 


dr. 

+  24-  +  *>J  = 


qj  1  dry  \ 


(2) 


where  Er^(r,f)  is  the  effective  local  electric  field  intensity 
and  B^r,f)  is  the  effective  local  magnetic  induction  field 
at  the  space-time  point  {r,f)  and  where  ry=r,(rsf)  de¬ 
scribes  the  displacement  of  the  electron  from  its  equilib¬ 
rium  position.  Here  qe  denotes  the  magnitude  of  the 
charge  and  m  the  mass  of  the  harmonically  bound  elec¬ 
tron  with  undamped  resonance  frequency  Wj  and  phenom¬ 
enological  damping  constant  $r  The  temporal  Fourier  in¬ 
tegral  representation  of  the  electric  and  magnetic  field 
vectors  of  the  effective  local  plane  electromagnetic  wave  is 
given  byfi 


(3) 
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when  magnetic  field  effects  are  neglected,  and 


/  2  2 

(*>  -  Uij 


{q3r/m2)E2ff 

2/ Sji i>)2  -  (qr!mc)'lE''ffu>1 


(18) 


is  defined  here  as  the  parallel  component  of  the  atomic  po¬ 
larizability,  where  ^,(o>t0)?0, 


3.  CRITICAL  FIELD  STRENGTH 

The  atomic  polarizability  is  then  seen  to  be  nonlinear  in 
the  local  electric  field  strength  when  magnetic  field  effects 
are  included.  The  nonlinear  term  appearing  in 

Eq.  (17)  will  be  negligible  in  comparison  with  the  linear 
term  when  the  Local  electric  field  strength  is  suffi¬ 

ciently  smaller  than  the  critical  field  strength  Ecrit 
=Ecnt(ti))  defined  by  the  condition 

(qJmcfE^ut2 

- z -  s  l, 

(<j2  -  u)j  +  2 i Sjtit)1  -  {qJm.c)*E2rlS<D2 

with  solution 

me  {ti>2  -  w2)2  +  A$ti>2 

Ecrit(tt>)  ^  A  &  2lV2  * 

\2qr  "  tf)  -  4  SjvfY1 

When  the  linear  term  in  Eq.  (17)  dominates 

the  nonlinear  term  so  that  ajX(wyEr^)  =*aj^(w).  Notice 
that  this  critical  electric  field  strength  depends  on  the 
value  of  the  applied  angular  frequency  to  of  the  electro¬ 
magnetic  wave  field.  Furthermore,  notice  that  this  critical 
field  strength  in  undefined  in  the  angular  frequency  band 
where 

w*  =  (w]  +  2<5f  ±  2  Sj  -Jwj  +  (20) 

for  each  resonance  feature  j,  where  wj  e  [wj\  w*].  Finally,  a 
rough  estimate  of  the  minimum  value  for  this  critical  field 
strength  is  seen  to  be  given  by 

me  hwf 

=  (21) 

fait 

for  each  resonance  feature,  where  fxB^qchi2mc  is  the 
Bohr  magneton.  The  quantity  hwj  is  recognized  as  the  en¬ 
ergy  level  state  of  the  harmonically  bound  electron, 

4*  NUMERICAL  EXAMPLE 

Consider  a  single  resonance  Lorentz  model  dielectric  (/ 
=0)  with  material  parameters  representative  of  the  near- 
infrared  line  in  triply  distilled  water  at  25  °Cy  where  te0 
=6.19  x  10 14  r/s  and  <5q-2.86  x  10 13  r/s*  The  angular  fre¬ 
quency  dependence  of  the  real  and  imaginary7  parts  of  the 
linear  atomic  polarizability  *  ajM  for  this  single  reso¬ 
nance  line  is  presented  in  Fig.  1,  and  the  angular  fre¬ 
quency  dependence  of  the  critical  field  strength  described 
by  Eq.  (19)  is  illustrated  in  Fig.  2,  where  ^=*5.91 
x  10u  r/s  and  wj  **  6.48  x  1014  r/s.  The  minimum  critical 
electric  field  strength  (Ecrit)mi„  =*  1.38  X  1011  V/m  is  seen 
to  occur  at  the  two  points  just  below  and  above  the  fre¬ 
quency  band  where  the  critical  field  strength  is 


undefined;  the  estimate  given  by  Eq.  (21)  yields  an  ap¬ 
proximate  value  of  7*5  x  10 11  V/m*  Notice  that  EcrU-+  «  at 
both  of  the  near-resonance  values  Wq  as  well  as  when  ei¬ 
ther  <i> — ►  0  or  cv— Notice  also  that  since  megavolt  per 
meter  electric  field  strengths  are  produced  by  ultrashort 
pulsed  ter  a  watt  laser  systems,  the  minimum  critical  field 
strength  for  this  example  is  at  least  five  orders  of  magni¬ 
tude  greater  than  that  which  is  currently  available. 

The  angular  frequency  dependence  of  the  relative  cor¬ 
rection  factor  for  the  perpendicular 

atomic  polarizability  given  in  Eq.  (17)  is  illustrated  in  Fig. 
3  for  the  near-infrared  resonance  line  of  triply  distilled 
water  for  several  values  of  the  relative  local  electric  field 
intensity  about  unity.  Notice  that  this  correction  factor 
reaches  its  peak  value  at  the  resonance  frequency  lo0 
when  \Eeff/(Eerit)mitt\<l  and  that  this  peak  value  bifur¬ 
cates  into  a  pair  of  symmetric  peaks  about  this  resonance 
frequency  when  \E€^i{Ecrit)min\  >  I,  a  local  minimum  now 


Fig,  1.  Real  (solid  curve)  and  imaginary  (dashed  curve)  parts  of 
the  linear  atomic  polarizability  for  the  near-infrared  resonance 
line  in  water 


Fig,  2.  Angular  frequency  dependence  of  the  critical  local  elec¬ 
tric  field  strength  (in  volts  per  meter)  for  the  near-infrared  reso¬ 
nance  line  in  water. 
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Fig,  6,  Angular  frequency  dependence  of  the  real  part  of  the  lon¬ 
gitudinal  atomic  polarizability  for  critical  and  subcritical  values 
of  the  local  electric  field  strength. 
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Fig.  7.  Angular  frequency  dependence  of  the  imaginary  part  of 
the  longitudinal  atomic  polarizability  for  critical  and  subcritical 
values  of  the  local  electric  field  strength. 


resonance  about  the  undamped  resonance  frequency  w0 
when  the  local  electric  field  strength  is  equal  to  or  less 
than  the  minimum  critical  field  strength  (EcrU)min  given 
in  £q,  (21),  However,  as  the  local  field  strength  exceeds 
this  minimum  critical  value,  this  single  resonance  peak  is 
found  to  split  into  a  pair  of  resonance  structures,  one 
downshifted  and  the  other  uphifted  in  angular  frequency 
from  the  characteristic  resonance  frequency  wo  of  the  Lor¬ 
entz  model  dielectric. 


5.  DISCUSSION 

The  results  presented  here  show  the  precise  manner  in 
which  the  magnetic  field  influences  the  atomic  polariz¬ 
ability  in  the  Lorentz  model  of  resonance  polarization.  In¬ 
clusion  of  the  magnetic  field  in  the  Lorentz  force  relation 
results  in  a  microscopic  polarization  density  that  contains 
both  perpendicular  and  parallel  components  relative  to 
the  plane  wave  propagation  vector  of  the  local  driving 
field  that  are  both  nonlinear  in  the  local  electrical  field 
strength.  This  nonlinearity  becomes  significant  when  the 
local  applied  electric  field  strength  exceeds  a  minimum 
critical  field  strength  whose  value  increases  linearly  with 
the  resonance  frequency  Numerical  calculations  show 
that  these  nonlinear  terms  are  entirely  negligible  for  ef¬ 
fective  field  strengths  that  are  typically  less  than 
^  10 12  V/m  and  that  they  begin  to  have  a  significant  con¬ 
tribution  for  field  strengths  that  are  typically  greater 
than  ^-1016  V/m  for  a  highly  absorptive  material.  Fortu¬ 
nately,  this  critical  field  strength  is  at  least  five  orders  of 
magnitude  greater  than  that  which  is  currently  available 
in  ultrashort  pulsed  terawatt  laser  systems.  Although 
this  magnetic  field  effect  is  negligible  in  most  practical 
situations,  it  has  been  shown  to  be  significant  in  materi¬ 
als  exhibiting  chirality,^11  and  it  may  also  become  impor¬ 
tant  for  artificial  materials.  Nevertheless,  seemingly  ex¬ 
cessive  electric  field  strengths  commonly  occur  in  nature; 
for  example,  it  is  estimated12  that  the  critical  field 
strength  for  the  alignment  of  water  molecules  for  crystal¬ 
lization  into  polar  ice  crystals  is  greater  than  IQ9  V/m. 


ftWj 


(24) 


hoij 

- — — . 
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(26) 


are  obtained* 

Similar  results  are  obtained  for  the  frequency  depen¬ 
dence  of  the  nonlinear  character  of  the  parallel  atomic  po¬ 
larizability  given  in  Eq.  (18).  Although  this  parallel  com¬ 
ponent  of  the  atomic  polarizability  vanishes  as  Eefj — >0,  it 
can  exceed  the  strength  of  the  transverse  component  of 
the  atomic  polarizability  when  Eeff  exceeds  the  minimum 
critical  field  strength.  The  frequency  dependence  of  the 
real  part  of  the  longitudinal  (parallel)  component  is  pre¬ 
sented  in  Fig*  6  and  the  imaginary  part  in  Fig*  7  for  two 
values  of  the  local  electric  field  strength.  Both  the  real 
and  the  imaginary  parts  are  seen  to  exhibit  a  pronounced 


Fig.  8.  Angular  frequency  dependence  of  the  relative  phasor  ve¬ 
locity  magnitude  of  the  harmonically  bound  electron  in  the  Lor¬ 
entz  model. 
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Abstract.  A  canonical  problem  of  central  importance  in  the  theory  of  ultrawideband  pulse  propaga¬ 
tion  through  temporally  dispersive,  absorptive  materials  is  the  propagation  of  a  Heaviside 
step-function  signal  through  a  medium  that  exhibits  anomalous  dispersion.  This  prob¬ 
lem  is  rich  in  the  use  of  asymptotic  theory,  SommerfaJd  and  Brillouin  provided  the  first 
(qualitatively  accurate  but  quantitatively  inaccurate)  doeed-form  approximations  of  the 
dynamic  evolution  of  this  waveform  through  a  single-resonance  Lorentz  model  dielectric 
based  upon  Debye's  method  of  steepest  descent.  An  improved  approximation  has  since 
been  provided  by  Oughstun  and  Sherman  using  modem,  uniform  asymptotic  methods 
that  rely  upon  the  aaddle-point  method.  An  accurate,  uniform  asymptotic  approxima¬ 
tion  describing  the  dynamical  evolution  of  the  unit  etep-f unction  modulated  aine  wave 
signal  through  a  single-resonance  Lorentz  model  dielectric  is  presented  here  baaed  upon 
their  work.  This  refined  asymptotic  description  results  in  a  continuous  evolution  of  the 
propagated  field  for  all  space-time  points. 

Kay  words,  asymptotic  methods,  dispersive  attenuative  wavs  propagation 

AMS  subject  classifications.  78M35,  78A40,  30E1B 
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I.  Introduction.  Asymptotic  analysis  is  widely  used  in  the  study  of  pulse  prop¬ 
agation  in  electromagnetics,  optics,  and  acoustics.  Typical  applications  include  the 
study  of  the  propagated  pulse  evolution  as  either  the  propagation  distance  or  the 
wavelength  becomes  either  large  or  small.  Here,  we  study  two-dimensional  (space  and 
time)  electromagnetic  pulse  propagation  through  an  unbounded  dielectric  material 
and  show  that  many  facets  of  uniform  asymptotic  theory  are  drawn  upon  in  order  to 
provide  a  continuous  asymptotic  approximation  to  the  propagated  pulse. 

Consider  a  linearly  polarized  plane- wave  electromagnetic  pulse  traveling  in  the 
positive  z-diiection.  The  material  through  which  the  pulse  travels  is  a  linear  di¬ 
electric  whose  relative  magnetic  permeability  ^  is  unity  and  whose  relative  dielectric 
permittivity  c(tu)  is  given  by  the  single-resonance  Lorentz  [IS]  dielectric  model  so  that 
the  complex  index  of  refraction  n(w)  —  is  given  by 

/  a  \  l/a 

(1.1)  n(w)  =  1  -  - £  —  )  . 

^  u2  —  ujg  +  2i6u  J 
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Brillouin  found  two  sets  of  saddle  points:  the  contributions  from  the  distant  saddle 
points  are  responsible  for  the  first  forerunner,  which  had  been  evaluated  for  «mail 
S  >  1  by  Sommerfeld;  the  contributions  from  the  near  saddle  points  give  rise  to  a 
second  forerunner  which  begins  after  the  arrival  of  the  first  forerunner.  Brillouin  also 
identified  the  main  signal,  character i2ed  by  its  harmonic  oscillation  at  the  carrier  fre¬ 
quency  we  of  the  input  pulse,  which  arrives  after  the  second  forerunner.  Sommerfeld 
and  Brillouin  mistakenly  concluded  that  both  forerunners  are  of  negligible  amplitude 
in  comparison  to  the  amplitude  of  the  main  signal. 

Although  Sommerfeld  and  Brillouin  correctly  predicted  the  existence  of  the  first 
and  second  forerunners,  they  did  not  accurately  describe  their  behavior.  It  wasn't  until 
1975  when  Oughstun  and  Sherman  [25,  26]  examined  Sommerfeld ’s  and  Brillouin *s 
seminal  works  that  an  accurate  description  of  these  two  forerunners  {now  known  as 
the  Sommerfeld  and  Brillouin  precursors)  and  the  main  signal  was  given  based  upon 
uniform  asymptotic  expansion  techniques.  Of  significant  importance,  Oughstun  [22] 
proved  that  the  Brillouin  precursor  can  be  of  significant  amplitude  when  the  rise  time 
of  the  pulse  is  faster  than  the  relaxation  time  of  the  material.  This  analysis  showed 
that,  asymptotically,  the  BriHouin  precursor  achieves  a  peak  amplitude  at  the  space* 
time  point  ct/z  —  ft(G),  where  n(0)  is  the  static  refractive  index  of  the  material,  and 
that  the  amplitude  of  this  point  decays  only  algebraically  as  z~ with  propagation 
distance  zy  whereas  the  amplitudes  of  the  Sommerfeld  precursor  and  the  main  mpml 
decay  exponentially*  The  relatively  alow  decay  rate  of  the  peak  amplitude  point  of 
the  Brhlouin  precursor  has  important  practical  applications  to  biomedical  imaging  as 
well  as  to  ground-  and  foliage-penetrating  radar  systems  (see  the  proceedings  series 
titled  Ultra-  Wideband  Shori-Polse  Electromagnetics  [1,  2,  7,  19,  34,  35])* 

The  exponential  appearing  in  (1*2)  may  be  written  in  terms  of  the  complex  phase 
function  defined  in  (1*3)  so  that,  with  some  algebraic  manipulation,  the  integral  rep¬ 
resentation  of  the  propagated  pulse  may  be  expressed  as  [29] 

(1.4)  exp  [f*(w.tf)] 

for  the  input  step-function  modulated  sine  wave  pulse  with  fixed  carrier  frequency 
>  0.  Because  of  the  exponential  appearing  in  the  integrand  of  (1*4),  the  Fourier- 
Lapiace  representation  is  ideally  suited  for  analysis  by  asymptotic  expansion  tech¬ 
niques  that  are  vakd  as  z  — t  oo  {i*e*,  as  the  propagation  distance  becomes  large), 
as  was  originally  done  by  Brillouin.  In  fact,  this  problem  is  much  richer  in  asymp¬ 
totic  theory  than  the  direct  application  of  the  method  of  steepest  descent.  In  order 
to  find  a  continuous  asymptotic  approximation  to  the  propagated  field  E(zyt),  three 
different  uniform  asymptotic  theories  are  required:  a  method  for  two  saddle  points 
symmetrically  located  about  the  imaginary  sods  whose  real  parts  are  located  at  plus 
and  minus  infinity;  a  method  for  two  first-order  saddle  points  that  coalesce  into  one 
second-order  saddle  point  at  some  fixed  space-time  point  and  then  separate  into  first- 
order  saddle  points  again;  and  finally,  a  method  for  first-order  saddle  points  in  the 
vicinity  of  a  simple  pole*  These  methods  have  been  used  extensively  by  Oughstun 
and  Sherman  to  study  pulse  propagation  through  Lorentz  model  dielectric  materials 
[21,  25  ,  26,  27,  28,  29];  their  results  form  the  basis  of  this  paper* 

The  material  parameters  (wo  =  4.0x1O1°/bi  =  20*0 xlO32/*3,  <5  =  0*28xlO1°/e) 

used  here  are  the  same  as  those  used  by  Sommerfeld  and  Brillouin  [6].  The  real 
and  imaginary  parts  of  the  complex  Index  of  refraction  n{u)  for  the  single-resonance 
Lorentz  model  dielectric  with  these  parameters  are  illustrated  in  Figure  1,1*  This 
choice  of  medium  parameters  corresponds  to  an  extremely  absorptive  medium.  Nev- 
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a  maximum  at  and -this  maximum  becomes  more  pronounced  ’as  z  oo*  Hence  *  -  - 
the  integral  appearing  in  (2*1)  may  be  approximated  by  the  contribution  due  to  a  - 
neighborhood  of  ,  the  accuracy  of  this  approximation  increasing  in  the  sense  of 
Poincare  [31]  as  z  — >  oo,  This  then  forms  the  principle  upon  which  the  saddle-point 
method  is  based.  The  saddle-point  method  may  be  directly  extended  to  integrals  of 
the  type  (see  Chapter  5  of  [29]) 

(2.2)  I{zt  0)  =  J  /( u>)  exp  [zgluf,  0)]  dw, 

in  which  the  position  of  the_saddla  point  of  the  exponential  function  g(w,  B)  is  a 
function  of  the  real- valued  parameter  0,  as  is  the  case  in  thus  paper*  This  is  achieved 
by  choosing  a  path  of  integration  which  moves  in  a  continuous  manner  as  6  varies 
continuously  over  some  specified  space-time  domain. 

3*  Behavior  of  the  Complex  Phase  Function.  In  order  to  apply  the  saddle-point' 
method  to  the  integral  appearing  in  the  integral  representation  (1,4)  of  the  propagated 
electric  field  component ’orast^p-fundtiM  modultLt*dsln¥  WEVe  'Sign&l,  the  behavior 
of  the  complex  phase  function  $(wp  0)  =  tw[Ti{tj)  ^~S\  musTT^Tcuo^Th^ lie  complex 
u-plane  [6T  29T  36]  as  a  function  of  the  real  parameter  0*  For  the  single-resonance 
Lorentz  model  dielectric  [18],  the  complex  phase  function  is  analytic  in  the 

iiz-plane  formed  by  the  two  branch  cuts  in  the  lower  half  of  the  u*-plane  symmetrically 
located  about  the  imaginary  axis.  Id  the  right  half  plane,  the  branch  cut  extends 
from  =  Vwq  -  -  1 5  to  s  -f-  w*  -  S2  -  i<5. 

The  saddle  points  of  ^  are  solutions  of  the  equation  ^(uip0)  =  0,  which  need 
only  be  solved  for  B  >  1  since  the  field  given  by  (1.4)  identically  vanishes  for  all 
0  Brillouin  [6]- first  showed  that  <f>  possesses -two  sets  of  saddle  points  and  pro¬ 
vided  first-order  approximations  of  the  locations  of  these  saddle  points;  more  accurate 
approximations  have  since  been  provided  by  Oughstun  and  Sherman  [29],  The  first 
set,  referred  to  as  the  distant  saddle  points  and  denoted  by  wsp±(0)»  consists  of  two 
first-order  saddle  points*  As  9  -t  1+,  ° 

(3-1)  lim  ujsp jt  cs  ±oo  -  i2 6, 

s_n+  SFu  1 

while  in  the  opposite  limit  as  0  — ►  oo,  the  distant  saddle  points  approach  the  outer 
branch  points 

(3.2)  \im^Msp±  S3  ±yju>%  4-  —  62'  —  %&  — 

respectively,  eo  that  |^SjP±(0)|  >  \/u^  +  u t*  for  all  0  >  1, 

The  second  set  of  saddle  points,  referred  to  as  the  near  saddle  points  and  denoted 
by  KspjJ  conslHtfi  of  two  first-order  saddle  points  that-  lie  on  the  imaginary  axis 
for  values  of  0  <  0Lr  coalesce  into  a  single  second-order  saddle  point  when  0  =  0 lf  and 
separate  into  two  first-order  saddle  points  symmetrically  located  about  the  imaginary 
axis  for  0  >  0It  where  01  «  0O  +  eo  that  9X  w  1,502  for  the  material 

parameters  used  here.  As  0  — ►  oo,  the  locations  of  the  near  saddle  points  approach 
the  inner  branch  points 

(3.3)  lim  a?  ±-s/u>g  —  S2  —  iS  ^ 

respectively.  Of  particular  interest  in  this  evolution  is  the  space-time  point  6  = 

60  =  n(0),  where  So  =  1.6  for  the  material  parameters  considered  here,  at  which  the 
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Fifi,  3.2  Contours  of  for  $  =  1,25  in  the  complex  u* -plane  along  with  an  acceptable 

deformed  path  of  intepnotixm  V(wi  Q),  Notice  that  the  original  and  deformed  contours  He 
on  the  same  side  of  the  pole * 


Fig.  3  J  Contours  of  0)}  for  $  m  9X  in  the  complex  v-plane  along  iwth  an  acceptable  de¬ 
formed  path  of  initiation  Notice  that  the  original  and  deformed  contours  lie  on 

the  same  side  of  the  pole. 
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introduced  a  change  of  variable  such  that  the  saddle  points  of  the  transformed  phase 
function  retain  the  limiting  values  of  the  original  saddle  points  and  then  expanded  the 
amplitude  function  as  a  finite  number  of  terms  plus  a  remainder  term,  the  remainder 
term  being  regular  and  equal  to  zero  at  the  saddle  points  of  the  transformed  phase 
function.  Here,  the  transformation  that  retains  the  essential  behavior  of  the  distant 
saddle  points  is  given  by 


(4.1) 


4 3 


where  and  0d($)  functions  to  be  determined.  The  amplitude  function  is 

then  written  as 


(4.2) 


•BoM), 


where  Hq{s,0)  is  regular  and  equal  to  zero  at  s  =  (±l/oy).  Substitution  of  (4.1) 
and  (4.2)  into  (L4)  results  in  the  uniform  asymptotic  description  of  the  Sommerfeld 
precursor  £s(zt£)  for  the  step-function  modulated  sine  wave  J29], 


Es(*,t)  ~  8  {exp  {-i-c0i{6)\  [70m  (|t*(9>) 

(4.3)  +2ati(0)e-i*/a71m  (f<*W))]}  -MM) 

asjr-^ocforalltf^l,  where  Jn{£)  denotes  the  Bessel  function  of  the  first  kind  of 
integer  order  n*  The  magnitude  of  the  remainder  term  is  bounded  as 

(4.4)  l*M)|  <  [|j,  (ad(8)^)  |  +  |j3  (cM<g)|)j] 

for  z  >  Z  >  0  and  6  >  1,  where  K  >  0  is  a  constant  independent  of  8  and  z.  The 
coefficients  appearing  in  (4.3)  are  given  by 


(4.5a) 

(4.5b) 


(4.5c) 


(4.5d) 


Ms)  =  \  [0(wSP+,9)  -  4>{uSP- , 9)]  =  -a  |d>(^SP+ , 9)}  , 
M«)  -  5  [tf(wSP+,9)  +  <K*sp-d , 0)}  =  at {^SP+J)} , 

;P+  -we  {.2ad(8))  (  i 


■»<*>  -  2 


*»iW 

(wSp+ 


nW  = 


a). 

_ i _ ( _ ( _ ) ,/! 

USP£  -  We  V  2o d{0)J  ^  i4>W(u>sp-,8)  J 

_ i _ ( 1 

U$P+  -we  \2ad{8)J  \i0<a)(wsp+,9) J 


i/a 


i  ( _ L-\  (  4ag(g)  A 

wsp-  “We  V  2ad{8) J  (  j^(2) (wSp. , 8) j 


i/s' 


For  values  of  0  bounded  away  from  1  such  that  (z/c)  |ad(0}[  >  1,  the  large  argument 
asymptotic  expansion  of  the  Bessel  function  may  be  substituted  into  (4.3).  This 
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would  be  bounded  by 

(5.1)  | An |  <  A\z\-V"+kV2  +  Beep  [R{*(*(Wj-,  8)  -  9(wu  5)]}) , 

where  Ay  B  are  some  constants  independent  of  zf  and  where  is  the  dominant  saddle 
point  while  Wj  is  the  other  saddle  point,  Le,,  8E{w;}  <  Because  u>4  is  the 

dominant  saddle  point,  the  second  term  in  (5,1)  is  negligible  in  comparison  to  the 
first,  for  large  enough  \z\.  However,  as  9  -+  0*,  ft{w4}  3R{t^}  and  \z\  must  increase 
without  bound  in  order  for  the  second  term  in  (5,1)  to  remain  negligible  compared  to 
the  first. 

A  uniform  expansion  of  the  Brillouin  precursor,  valid  for  ah  0  >  1,  is  obtained 
through  use  of  the  theorem  [9,  12,  29]  originally  due  to  Chester,  FViedman,  and 
UrseU  [9],  Here,  the  transformation  that  retains  the  behavior  of  the  saddle  points 
about  $  =  Bi  is  given  by 

(5.2)  j*'3  -  ai{&)v  ~  ao(9)  +  =  0, 

where  Oo,ai  are  functions  to  be  determined.  Note  that  there  are  three  possible 
branches  of  the  inverse  function.  Chester,  Friedman,  and  Ursell  [9]  proved  that  only 
one  branch  of  the  transformation  defines  a  conformal  mapping  of  some  disc  that 
contains  both  saddle  points.  It  is  this  branch  which  must  be  used,  as  shown  by 
Bleistein  and  Handelsman  [4].  As  before,  the  amplitude  function  is  expanded  as  a 
finite  number  of  terms  plus  a  remainder  that  is  regular  and  equal  to  zero  at  the 
transformed  saddle  points.  This  is  obtained  using  the  expansion 

(5.3)  — —  £  s  G0(M)  =  MO)  +  M*)"  +  O'3  -  ai(*))Ao(M), 

—  U)c  (XV 

where  and  Hq  are  to  be  determined.  Substitution  of  (5.2)  and  (5,3)  into  (1.4) 

results  in  the  uniform  asymptotic  expansion  of  the  Brihouin  precursor  £a(z,t)  for 
the  step- function  modulated  sine  wave  (29] 


EB(z,t)  =  -R  j  exp  [~oo(e]  |  (^)V3  e_i3w/3  ■  Ai  (ttl{ff)  e~t2''s  (|)V8) 

(i  [ - - — -M*)  + - - - fc-(0)l  +o(i)} 

l 2  rsp+  -  «*  ^sp-  “w<=  J  w  J 

+  (£)2/3  e-^/3^0)  «-*-/»  Qa/S) 


(5.4) 


{*>!*( 


(0)  [U>SP+  ~  We 


h+(8)~ 


WSPZ  ~  u><- 


M*>  + 


-(;)}}}• 


where  Ai(()  denotes  the  Airy  function  and  the  coefficients  appearing  in  (5.4)  are  given 
by 


(5,5a) 


£*o{0)  =  ~  [l)KwSP+ .  0)  +  ‘t’&SP-  *  ^)]  ■ 


(5.5b) 

(5.5c) 


M*)  = 


/  2a\'\0)  \ 
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total  field  evolution  as  z  — t  oo.  The  instantaneous  angular  frequency  of  oscillation 
of  the  Brillouin  precursor  starts  at  zero  when  0  —  0i  and  monotonically  increases  to 
approach  the  value  y/tJ^  —  S2  with  increasing  0. 

Previous  work  [29,  36)  exhibited  a  discontinuity  in  the  uniform  asymptotic  ap¬ 
proximation  of  the  Brillouin  precursor  about  the  space- time  point  9  —  0* ,  which  the 
authors  attributed  to  numerical  instabilities.  These  numerical  instabilities  are  now 
known  to  be  caused  by  the  use  of  an  unnecessary  approximation  of  the  lower  near 
sad  die- point  location  ^SP-  {9)  and  incorrect  phase  values  for  the  coefficients  kt(0)  for 
all  9  >  $i-  Here,  all  saddle-point  locations  and  coefficients  are  numerically  determined 
so  that  this  approximation  error  is  avoided,  resulting  in  a  continuous  evolution  of  the 
Brillouin  precursor  for  all  6  >  1. 

6.  The  Signal  Contribution,  The  signal  contribution  to  the  propagated  field  of 
the  step-function  modulated  sine  wave  is  due  to  the  simple  pole  singularity  at  ui  ~  w,. 
appearing  in  the  integrand  of  (1.4).  It  is  assumed  here  that  uc  is  real,  positive, 
and  finite  (as  all  physical  frequencies  must  be).  The  simple  pole  located  at  may 
influence  the  value  of  E{z%  £)  if  the  path  ?(w,  9)  crosses  the  pole  singularity,  or  if  either 
of  the  saddle  points  uijp+  ^  (0)  comes  within  close  proximity  of  the  pole,  or  both. 

Let  the  path  V(uj}$)  comprise  the  portions  of  the  steepest  descent  paths  em¬ 
anating  from  and  and  continuing  into  the  upper  half  plane.  In 

the  deformation  of  the  original  path  of  integration  appearing  in  (L4)  Into  the  path 
V(u;)  0),  the  simple  pole  located  at  wc  is  crossed  at  some  space- time  point  0  =  9t  >  1* 
Hence,  the  original  integral  and  the  integral 


(6.1) 

EsDp(z,t)  =  „  s  < »  /  exp 

2*  {  /p(w,9)W-Wc 

[^(w.fl)]  dwj 

along  the  deformed  contour  are  related  by 

(6.2a) 

E(z,t)  =  Esdp{z  ,t) 

for  0  <  0#1 

(6.2b) 

E(z,t)  =  ESDp(z,t)  - 

for  0  =  9i3 

(6.2c) 

E(z,t)  —  E$dp{zi  0  —  2ir*7e^e^“‘”^ 

for  0  >  9m% 

where 

(6.3) 

7  =  lim  (w  —  wc) — - —  = 

U— U)  — 

■  t 

is  the  residue  of  the  amplitude  function  at  the  simple  pole  u?c.  Direct  application  of 
the  saddle-point  method  to  yields  a  uniform  expansion.  However,  it  does 

not  provide  a  uniform  expansion  for  J2(2,£),  In  order  to  see  this,  let  0C  denote  the 
space-time  point  at  which  the  pole  becomes  dominant  over  the  saddle  point 
in  the  sense  that  0)}  >  &{^{a>sF,  0)}.  Because  the  path  ^(^,0)  lies  in 

the  valleys  of  u>Spt  it  necessarily  follows  that  0#  <  0C>  Prom  (6.2),  it  is  obvious 
that  E(z3t)  changes  discontinuously  about  0  =  0*  when  the  saddle- point  method  is 
applied  to  Esnp{zjt)*  However,  because  the  simple  pole  lies  in  the  valley  of  the 
saddle  point,  the  residue  contributions  appearing  in  (6.2)  are  exponentially  smaller 
than  Esdp(z » 0  &nd  pose  no  problem  in  obtaining  a  uniform  approximation  to  E(z ,  f). 
This  is  not  the  case  at  0  =  9C  when  the  simple  pole  wc  becomes  dominant  over 
the  saddle  point  usp-  Here,  &{d>(wsp,  0)}  <  0)}  for  0,  <  0  <  0C  and  the 
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while  &§  and  the  appropriate  residue  .contribution  are  determined  by, the -saddle  point 
whose  path  of  steepest  descent  crosses  the  simple  pole  located  ,  at  .The  space- time 
point  &M  at  which  the  path  0)  crosses  the  simple  pole  located  at  is  determined 
by  the  equation  *  - 

(6-7) 


where  Y  m  and  u$p  denotes  either  or  wsp+(0).  We  have  found  (by 

inspection)  that  the  deciding  factor  as  to  which  saddle  point  determines  the  value 
of  B9  is  the  value  of  Y( wCT  0}  at  $  =  1.  If  Y  (wc,  1)  >  G,  the  steepest  descent  path 
emanating  from  the  near  saddle  point  wsp+(0)  is  used  to  determine  the  value 
If  y(u>CJ  1)  <  0,  the  steepest  descent  path  emanating  from  the  distant  saddle  point 
w£P+  (#)  is  used  to  determine  the  value  Bt,  If  Y(u)Cf  1)  s=  0,  then  B9  =  L  The  value  of 

wc  =  (tty  which  separates  the  two  cases  is  given  approximately  by  u ty  **  4*2925  x  1016 
rad/s  for  the  material  parameters  considered  here.  As  a  consequence,  we  have  applied 
the  uniform  theoiy  to  account  for  two  first-order  ©addle  points  with  a  nearby  simple 
pole  singularity*  As  an  example,  for  angular  frequencies  satisfying  Y(u)€l  1)  >  0, 
the  uniform  signal  contribution  is  given  by 


=  d-Si  ji7tl  |-iir  erfc  (iAD(0)y^  exp  [|*{«e,0) J 

+ZFwV^“p[^(“s,,S'S)]] 

+  »7v  -iff  erfc  exp  ^{wci  0)] 

<6>8a>  +A^(?)VTexp[^{w^'?)]]}'  6<$ - 

Ee(2,t)  =  d-ffi  j*7u  ^  iff  erfc  (iAD(fi)y^  exp  [^{we,fl)J 

+  t7„  iff  erfc  (-iA* exp  |~4>(wei0)J 

{6,8b)  +Ajv{^)^?eXP  }  ’i"®  {*exP  [~^(we.  0)]}  > 

as  z  — ►  oo*  Here,  7^  =  i  is  the  residue  of  the  amplitude  function  i(w  —  wc)  at 
erfc(£)  is  the  complementary  error  function,  and 


wc* 


(6,9)  Aj p>N(0)s  \<f>{usp+'„  1 ~ ^)]  '  • 

A  similar  expression  may  be  given  for  the  uniform  signal  contribution  for  carrier 
frequencies  wc  satisfying  Y(u>Ci  1)  >  0  [8],  _ 

Questions  about  the  accuracy  of  this  uniform  expansion  as  it  applies  to  the  clas¬ 
sical  problem  of  a  step-function  modulated  sine  wave  signal  have  persisted  because, 
until  now,  it  has  been  difficult  to  isolate  the  leading  edge  of  the  pole  contribution 
from  the  remainder  of  the  field*  For  all  values  of  either  the  Sommerfeld  or  the 
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Fig,  6.2  Comparison  of  asymptotic  (dashed  curve)  and  numerical  (solid  curve)  pole  contributions 
for  a  step-function  modulated  sine  want  unt4  intro- absorption  band  frequency  usc  =  1.25cjq 
at  an  observation  distance  of  s  ra  5.3*^.  The  asymptotic  apprtHarftation  is  made  (a)  uti¬ 
lizing  only  the  distant  saddle  point  w,jPt{0),  and  (b)  utilising  both  saddle  points  uj_p+  (0) 

andu3P+W'  5  D 


at  a  propagation  distance  of  z  &  5.3 zd.  In  Figure  6.2  (b),  the  asymptotic  approxi¬ 
mation  of  the  pole  contribution  (dashed  curve)  is  calculated  using  both  msp+  (0)  and 
The  open  circle  denotes  the  space-time  value  $f  when  the  path 
crosses  wc.  A  comparison  of  the  two  figures  shows  that  the  contribution  from  the  near 
saddle  point  msp+  (0)  is  minimal.  However,  these  figures  again  confirm  that  the  uni¬ 
form  asymptotic  theory  provides  the  correct  approximation  to  the  pole  contribution. 
This  is  in  contrast  to  the  previous  hypothesis  [29]  that  both  saddle  points  u;5i3+(0) 
and  (0)  are  sufficiently  removed  from  uc  when  wc  is  within  the  absorption  band 
of  the  material  so  that  i?c(2,  £}  is  composed  solely  of  the  residue  term  which  abruptly 
begins  at  6  —  BM. 

7.  Conclusions.  The  propagated  electric  field  component  of  a  linearly  polarized, 
plane- wave,  step-function  modulated  sine  wave  signal  of  applied  carrier  frequency 
wc  traveling  through  a  single- resonance  Lorentz  dielectric  has  been  studied  in  this 
paper.  An  analytic,  asymptotic  approximation  of  the  propagated  wave  field  requires 
the  application  of  three  uniform  theories  in  the  subject  of  asymptotic  expansions 
of  integrals.  The  results  presented  here  provide  the  first  correct  description  of  the 
pole  contribution  of  the  field  (1.4)  for  applied  carrier  frequencies  that  lie  within  and 
above  the  absorption  band  of  the  material.  Previous  attempts  to  isolate  the  pole 
contribution  from  the  rest  of  the  field  were  unsuccessful  due  to  inaccurate  asymptotic 
results  for  the  Brillouin  precursor  about  the  space-time  point  0  *  0i  and  an  incorrect 
asymptotic  representation  of  the  pole  contribution  for  applied  carrier  frequencies  in 
the  range  u/c  >  Vwjj  +  ut*  —  P. 

Four  different  carrier  frequencies  of  the  step-function  modulated  sine  wave  signal 
were  considered  in  this  paper:  the  beiow-resonanee  carrier  frequency  uc  =  0.5^o,  the 
intra-absorption  band  carrier  frequencies  ujc  —  u>o  and  wc  =  !-25wo,  and  the  above- 
resonance  carrier  frequency  wc  —  2.5u>0.  The  total  asymptotic  field,  which  is  the 
sum  of  the  three  asymptotic  components,  £?(z,f)  =  +  7?n(-z,t)  +  Bc{z » *)t 

is  presented  for  each  applied  carrier  frequency  as  the  dashed  curves  in  the  upper 
plots  of  Figures  7.1-7.4.  The  solid  curves  appearing  in  these  figures  represent  the 
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e  =  ct/z 


Pig,  7J  TAe  total  asymptotic  (dashed  curve)  and  ntitnericaf  (solid  curve)  fields  of  the  step- 
function  modulated  me  wave  signal  with  intro- absorption  band  applied  carrier  frequency 
Wj  =  I>35u*d  rad/s  at  a  distance  of  x  =  3  X  10"*  m  »  5.3*^  into  the  single-resonance 
Lorentx  dielectric.  The  lower  figure  shows  the  difference  between  tAe  two  reeufcj. 


®  =  Q t/Z 


Fig.  7.4  TAe  total  asymptotic  (dashed  curve)  and  nnmericai  (solid  curve)  fields  of  the  step- function 
modulated  sine  waue  ei^naf  with  above- absorption  band  applied  carrier  frequency  u/c  —  2.&wq 
at  a  dutance  of  x  =  Q  X  I0”r  m  at  2.7*^  into  tAe  ftn^ie-rejonance  Lorcnte  dielectric .  The 
lower  figure  shows  the  difference  between  tAe  two  results 
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